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Abstract 

A cosmologically significant population of very luminous high-redshift galaxies 
has recently been discovered at submillimeter (submm) wavelengths. Advances in 
submm detector technologies have opened this new window on the distant Universe. 
Here we discuss the properties of the high-redshift submm galaxies, their significance 
for our understanding of the process of galaxy formation, and the selection effects 
that apply to deep submm surveys. The submm galaxies generate a significant frac- 
tion of the energy output of all the galaxies in the early Universe. We emphasize 
the importance of studying a complete sample of submm galaxies, and stress that 
because they are typically very faint in other wavebands, these follow-up observa- 
tions are very challenging. Finally, we discuss the surveys that will be made using 
the next generation of submm- wave instruments under development. 
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1 Introduction 



Discovering the process by which the dense, gravitationally bound galaxies 
formed in the Universe from an initially almost uniform gas, and understand- 
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ing the way their constituent populations of stars were born is a key goal of 
modern physical cosmology. A wide range of well understood physical pro- 
cesses are involved; including general relativity, gas dynamics and cooling 
physics, nuclear reactions and radiative transfer. However, the range of pos- 
sible initial conditions and the non-linear nature of most of the events, start- 
ing with the collapse of primordial density perturbations, ensure that these 
intimately connected processes can generate a very wide range of possible 
scenarios and outcomes. Galaxy formation can be studied by attempting to 
reproduce the observed Universe via analytical models and numerical simula- 
tions. The information required to constrain these models is provided by both 
forensic studies of the current constituents of the Universe, including stellar 
ages, chemical abundances and the sizes and shapes of galaxies, and by direct 
observations of the galaxy formation process taking place in the young Uni- 
verse at great distances. Direct observations exploit both the light emitted by 
distant galaxies, and the signature of absorption due to intervening structures 
along the line of sight, and began almost 50 years ago using sensitive optical 
and radio telescopes. Astronomers must now use all available frequencies of 
radiation to probe the properties of the Universe, from the lowest energy radio 
waves to the highest-energy 7-rays. It is vital to combine the complementary 
information that can be determined about the constituents of the Universe at 
different wavelengths in order to make progress in our understanding. 

This review discusses the results of a new type of direct observation of the 
galaxy formation process, made possible by the development of powerful new 
radiation detectors sensitive to wavelengths in the range 200 /im to about 
1mm: the submillimeter (submm) waveband. The detection of submm radi- 
ation from distant galaxies is one of the most recent developments in obser- 
vational cosmology, and has finally brought this region of the electromagnetic 
spectrum into use for making cosmological observations not directly connected 
with the cosmic microwave background (CMB; Partridge and Peebles, 1967). 
With the possible exception of the hardest X-ray wavebands, studies of distant 
galaxies in the submm waveband remained elusive for the longest period. We 
will also discuss some observations at the mid- and far-infrared(IR) wavebands 
that bound the submm waveband at short wavelengths, usually defined as the 
wavelength ranges from about 5-40 and 40-200 /im, respectively. 

The most significant reason for the late flowering of submm cosmology is the 
technical challenge of building sensitive receivers that work efficiently at the 
boundary between radio-type coherent and optical-like incoherent detection 
techniques. In addition, atmospheric emission and absorption permits sensi- 
tive submm observations from only high mountain sites, and only in specific 
atmospheric windows. The zenith opacity from the best sites in the clearest 
submm atmospheric window at 850 /im is typically about 0.1. Furthermore, 
the long wavelength of submm radiation limits spatial resolution unless very 
large filled or synthetic apertures are available. The largest single apertures 
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available at present are in the 10-30 m class, providing spatial resolution of 
order 10 arcsec. This resolution is much coarser than the sub-arcsec resolu- 
tion of optical and near-IR observations. The appearance of the same region 
of sky at optical and submm wavelengths is compared in Fig. 1 to illustrate 
this point: the multicolor optical image was obtained using the Hale 5-m tele- 
scope at Mt. Palomar, while the 850-/zm submm image was obtained using 
the 15-m James Clerk Maxwell Telescope (JCMT) on Mauna Kea. Interfer- 
ometers can dramatically enhance the resolution of images, but so far have 
only operated at longer mm wavelengths. The commissioning of the 8-element 
Sub-Millimeter Array (SMA; Ho, 2000)0 on Mauna Kea in Hawaii with base- 
lines of up to about 500 m, the first dedicated submm- wave interferometer, will 
provide images with subarcsecond resolution. The much larger 64-element At- 
acama Large Millimeter Array (ALMA; Wootten, 2001)0 will be in service at 
the end of decade. 



A key development was the commissioning of the Submillimetre Common- 
User Bolometer Array (SCUBA) camera at the JCMT in 1997 (Holland et al., 
1999). SCUBA images the sky in the atmospheric windows at both 450 and 
850 /iin in a 2.5-arcmin-wide field, using hexagonal close-packed arrays of 91 
and 37 bolometer detectors at the respective wavelengths. SCUBA provided 
a dramatic leap forward from the pre-existing single-pixel or one-dimensional 
array instruments available. The combination of field of view and sensitivity 
was sufficient to enable the first searches for submm-wave emission from previ- 
ously unknown distant galaxies. The Max-Planck Millimetre Bolometer Array 
(MAMBO; Kreysa et al., 1998) is a 1.25-mm camera with similar capabili- 
ties to SCUBA, which operates during the winter from the Institut de Radio 
Astronomie Millimetrique (IRAM) 30-m telescope on Pico Veleta in Spain. A 
similar device — the SEST Imaging Bolometer Array (SIMBA) — designed at 
Onsala in Sweden is soon to begin operation on the 15-m Swedish-ESO Sub- 
millimetre Telescope (SEST) in Chile, providing a sensitive submm imaging 
capability in the South. The capability of mm and submm-wave observatories 
is not standing still: a number of larger, more sensitive mm- and submm- 
wave cameras are under construction, including the SHARC-II (Dowell et al., 
2001), BOLOCAM (Glenn et al., 1998) and SCUBA-II instruments^ Bolome- 
ter technology continues to advance. The advent of extremely stable super- 
conducting bolometers that require no bias current and can be read out using 
multiplexed cold electronics, should ultimately allow the construction of very 
large submm detector arrays of order 10 4-5 elements (for example Benford et 
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4 Details can be found in Table 3. The next-generation SCUBA-II camera for the 
JCMT is under development at the United Kingdom Astronomy Technology Cen- 
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Fig. 1. A comparison of deep optical and submm views of the sky. The back- 
ground image is a 3-color optical image of the rich cluster of galaxies Abell 1835 
at the low/moderate redshift z = 0.25 (Smail et al., 1998b) taken using the 5-m 
Hale telescope, overlaid with the 14-arcsec resolution contours of a SCUBA 850-/xm 
submm- wave image of the same field (Ivison et al., 2000a). North is up and East to 
the left. The brightest SCUBA galaxies at (-45,-15), (65,0) and (20,-60), and the 
central cD galaxy (Edge et al., 1999), all have clear radio detections at a frequency 
of 1.4 GHz in images with higher spatial resolution than the SCUBA contours, ob- 
tained at the Very Large Array (VLA), supporting their reality. The bright SCUBA 
galaxy at (—45, —15) is associated with SMM J14011+0253, an interacting pair of 
galaxies at redshift z = 2.56 in the background of the cluster (Frayer et al., 1999). 
Spectacular fragmented structure appears in the Easterly red component of this 
galaxy in Hubble Space Telescope (HST) images (Fig. 18). 

al., 2000). SCUBA-II is likely to be the first instrument to exploit this tech- 
nology, providing a 8 x 8-arcmin 2 field of view at the resolution limit of the 
JCMT. 

The first extragalactic submm/mm surveys using SCUBA and MAMBO re- 
vealed a population of very luminous high-redshift galaxies, which as a pop- 
ulation, were responsible for the release of a significant fraction of the energy 
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generated by all galaxies over the history of the Universe (Blain et al., 1999b). 
Almost 200 of these galaxies are now known (Smail et al., 1997; Barger et al., 
1998; Hughes et al., 1998; Barger et al., 1999; Eales et al., 1999, 2000; Lilly 
et al., 1999; Bertoldi et al., 2000; Borys et al., 2002; Chapman et al., 2002a; 
Cowie et al., 2002; Dannerbauer et al., 2002; Fox et al., 2002; Scott et al, 
2002; Smail et al., 2002; Webb et al., 2002a). There is strong evidence that 
almost all of these galaxies are at redshifts greater than unity, and that the 
median redshift of the population is likely to be of order 2-3 (Smail et al., 
2000, 2002). However, only a handful of these objects have certain redshifts 
and well-determined properties at other wavelengths (Frayer et al., 1998, 1999; 
Ivison et al., 1998, 2001; Kneib et al., 2002). The results of these mm/submm 
surveys provide complementary information to deep surveys for galaxies made 
in the radio (Richards, 2000), far-IR (Puget et al., 1999), mid-IR (Elbaz et 
al., 1999) and optical (Steidel et al., 1999) wavebands. Submm observations 
are a vital component of the search for a coherent picture of the formation 
and evolution of galaxies, which draws on data from all wavebands where the 
distant Universe can be observed. 

In this review, we describe the key features of the submm emission processes in 
galaxies. We summarize the current, developing state of submm- wave observa- 
tions of distant galaxies, including the results of both blank-field surveys, and 
targeted observations of known high-redshift galaxies, including radio-galaxies, 
optically-selected quasars/QSOs, X-ray detected active galactic nuclei (AGNs) 
and optically-selected Lyman-break galaxies (LBGs). Submm- wave surveys 
are not immune to selection effects, and we discuss their strengths and weak- 
nesses. We describe the properties of the class of submm-luminous galaxies, 
and discuss the key results that are required to make significant progress in un- 
derstanding them. We consider the relationship between the submm-selected 
galaxies and other populations of high-redshift galaxies, and describe models 
that can account for the properties of submm-selected galaxies. We introduce 
the unusually significant effects of the magnification of distant submm-selected 
galaxies due to gravitational lensing (Schneider et al., 1992). Finally, we recap 
the key developments that are keenly awaited in the field, and describe some 
of the exciting science that will be possible in the next decade using future 
instruments. 

The cosmological parameter values assumed are generally listed where they 
appear. We usually adopt a flat world model with a Hubble constant H = 
65kms _1 Mpc -1 , a density parameter in matter fi m = 0.3 and a cosmological 
constant = 0.7. 
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2 Submm-wave emission from galaxies 

There are two major sources of submm radiation from galaxies: thermal con- 
tinuum emission from dust grains, the solid phase of the interstellar medium 
(ISM), and line emission from atomic and molecular transitions in the inter- 
stellar gas. The ladder of carbon monoxide (CO) rotational transitions, spaced 
every 115 GHz, is the most important source of molecular line emission, but 
there is a rich zoo of other emitting molecules in the denser phases of the 
ISM. Submm surveys for distant galaxies have so far been made using cam- 
eras that detect only continuum dust emission, and so this will be the main 
focus of this review. However, the search for line emission is already important, 
and its study will become increasingly significant. The spectral resolution pro- 
vided by line observations reveals much more about the physical and chemical 
conditions in the ISM, for studies of kinematics, metallicity and excitation 
conditions. Molecular lines can also be used to obtain a very accurate spectro- 
scopic redshift for the ISM in high-redshift galaxies with prior optical redshifts 
(for example Frayer et al., 1998). Searches for redshifts at cm and (sub)mm 
wavelengths using CO lines will be possible using future telescopes. 

The best studied regions of the Universe in the submm waveband are Giant 
Molecular Clouds (GMCs) in the Milky Way, in which ongoing star formation 
is taking place (Hollenbach and Tielens, 1997). GMCs are perhaps very low- 
luminosity archetypes for distant dusty galaxies, although these galaxies have 
far-IR luminosities that are up to 4 orders of magnitude greater than that of 
the whole Milky Way. 

Detailed, resolved submm-wave images and spectra only exist for low-redshift 
galaxies (for example Regan et al., 2001; Sakamoto et al., 1999), and it is often 
necessary to use them as templates to interpret the properties of more distant 
galaxies. A very important class of well-studied galaxies similar in luminosity, 
and perhaps in physical properties, to high-redshift submm galaxies are the 
ultraluminous IR galaxies (ULIRGs) discovered in the InfraRed Astronomy 
Satellite (IRAS) all-sky survey in the mid 1980's (see the review by Sanders 
and Mirabel, 1996). ULIRGs are usually defined as having a bolometric lumi- 
nosity, integrated over all wavelengths at which dust emission dominates the 
SED (from about lmm-8/im), in excess of 10 12 L Q F£] They are amongst the 
most luminous of all galaxies, but number less than 0.1% of galaxies in the local 
Universe. Due to their selection by IRAS, they are typically at relatively low 
redshifts, less than about 0.3. The first JiMS'-detected high-redshift ULIRG 
was identified by Rowan-Robinson et al. (1991) at z = 2.3. The current record 
redshift for a galaxy detected by IRAS is z = 3.9 for APM 08279+5255 (Irwin 
et al., 1998). Both these galaxies appear to be extremely luminous; however, 
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Fig. 2. Various observed restframe spectral energy distributions (SEDs) of galax- 
ies from the radio to the near-IR wavebands. Two examples of the most lumi- 
nous low-redshift galaxies detected by IRAS are included (I). Five very luminous 
high-redshift galaxies that have been, or could have been, detected directly in deep 
submm surveys (S), three high-redshift galaxies serendipitously magnified and made 
easier to study by the gravitational lensing effect of foreground galaxies and also 
detected by IRAS (L), and five high-redshift AGNs detected in optical or radio 
surveys (H) are also shown. In addition, three template SEDs are shown. One in- 
cludes the properties of CO and atomic fine-structure emission lines in the (sub)mm 
waveband at wavelengths from 100 to 3000 /im (Blain et al., 2000b), one includes 
polycyclic aromatic hydrocarbon (PAH) molecular emission features at wavelengths 
~ 10 /im in the mid-IR waveband (Guiderdoni et al., 1998), and one is normalized 
to the typical SED of a sample of low-redshift IRAS galaxies (Dunne et al., 2000). 
For further information on far-IR SEDs see Dale et al. (2001). With the exception 
of the high-redshift AGNs and the lensed galaxies, the templates tend to provide a 
reasonable description of the SED at wavelengths around and longer than its peak, 
the regime probed by submm surveys. Less luminous galaxies like the Milky Way 
have dust spectra that peak at a wavelength about a factor of 2 longer than these 
templates (Reach et al., 1995). 

their luminosities are boosted by at least a factor of ten due to gravitational 
lensing by foreground galaxies. A compilation of the properties of some of the 
most extreme ULIRGs is given by Rowan- Robinson (2000). The IR spectral 
energy distributions (SEDs) of some low-redshift ULIRGs and a compilation 
of results for the more sparsely sampled SEDs of high-redshift dusty galaxies 
are illustrated in Fig. 2. 
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2. 1 The power source for dusty galaxies 

About 99% of the energy released by galaxies in the submm and far-IR wave- 
bands is produced by thermal emission from dust grains; the remainder comes 
from fine-structure atomic and molecular rotational line emission. However, 
the source of the energy to power this emission by heating dust is often un- 
clear. Any intense source of optical/ultraviolet (UV) radiation, either young 
high-mass stars or an accretion disk surrounding an AGN, would heat dust 
grains. Because dust emits a featureless modified blackbody spectrum, submm 
continuum observations can reveal little information about the physical con- 
ditions within the source. Regions of intense dust emission are very optically 
thick, and so little information can be obtained by observing optical or UV 
radiation. 

In typical spiral galaxies, with relatively low far-IR luminosities of several 
1O 1O L (for example Alton et al., 2000, 2001), the dust emission is known to 
be significantly extended, on the same scale as the 10-kpc stellar disk.0 The 
emission is certainly associated with molecular gas rich star-forming regions 
distributed throughout the galaxy (Regan et al., 2001), in which dust is heated 
by the hot, young OB stars. 

In intermediate luminosity galaxies, such as the interacting pair of spiral galax- 
ies NGC 4038/4039 'the Antennae' (Mirabel et al., 1998; Wilson et al., 2000), 
the most intense knots of star-formation activity, from which most of the lu- 
minosity of the system emerges, are not coincident with either nucleus of the 
merging galaxies, but occur in a deeply dust-enshrouded overlap region of the 
ISM of the galaxies. This provides a strong argument that almost all of the 
energy in this system is being generated by star formation rather than an 
AGN. 

In more luminous ULIRGs that are at sufficiently low redshift for their internal 
structure to be resolved, the great majority of the dust emission arises in a 
much smaller, sub-kpc region (Downes and Solomon, 1998; Sakamoto et al., 
1999) within a merging system of galaxies. It is plausible that a significant 
fraction of the energy could be derived from an AGN surrounded by a very 
great column density of gas and dust that imposes many tens of magnitudes 
of extinction on the emission from the AGN in the optical and UV wavebands, 
and which remains optically thick even at near-IR wavelengths. Alternatively, 
an ongoing centrally condensed burst of star-formation activity, fueled by gas 
funneled into the center of the potential well of a pair of interacting galaxies 
by a bar instability (Mihos, 2000) is an equally plausible power source. 

If the geometry of absorbing and scattering material is known or assumed, 
6 1 pc = 3.09 x 10 16 m 
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then radiative transfer models can be used to predict the SED of a galaxy, 
which should differ depending on whether the source of heating is a very 
small AGN with a very hard UV SED, or a more extended, softer-spectrum 
nuclear star- forming region (for example Granato et al., 1996). Note that the 
results are expected to be very sensitive to the assumed geometry (Witt et 
al., 1992). In merging galaxies this geometry is highly unlikely to be spherical 
or cylindrical, and is uncertain for the high-redshift galaxies of interest here. 
In the case of AGN heating, the SED would be is expected to peak at shorter 
wavelengths and the mid-IR SED would be expected to be flatter as compared 
with a more extended star-formation power source. Both these features would 
correspond to a greater fraction of hot dust expected in AGNs (see Fig. 2), 
and is seen clearly in the SEDs of low-redshift IRAS-detected QSOs (Sanders 
and Mirabel, 1996). 

An alternative route to probing energy sources in these galaxies is provided by 
near- and mid-IR spectroscopy. At these longer wavelengths, the optical depth 
to the nucleus is less than in the optical/UV, and so the effects of the more 
intense, harder UV radiation field expected in the environs of an AGN can be 
observed directly. These include the excitation of characteristic highly-ionized 
lines, and the destruction of relatively fragile polycyclic aromatic hydrocarbon 
(PAH) molecules (Rigopoulou et al., 1999; Laurent et al., 2000; Tran et al., 
2001), leading to the suppression of their distinctive emission and absorption 
features. Mid-IR spectroscopic observations with the successor to IRAS, the 
Infrared Space Observatory (ISO) in the mid 1990 's indicated that most of the 
energy from low-redshift ULIRGs is likely generated by star-formation activity 
rather than AGN accretion. However, the fraction of ULIRGs containing AGN 
appears to increase at the highest luminosities (Sanders, 1999). This could be 
important at high redshifts, where the typical luminosity of dust-enshrouded 
galaxies is greater than in the local Universe. In addition, there may be duty- 
cycle effects present to make an AGN accrete, and perhaps to be visible, for 
only a fraction of the duration of a ULIRG phase in the evolution of the galaxy 
(Kormendy and Sanders, 1988; Sanders et al., 1988; Archibald et al., 2002). 
X-ray observations also offer a way to investigate the power source, as all but 
the densest, most gas-rich galaxies, with particle column densities greater than 
10 24 cm~ 2 are transparent to hard (> 2keV) Xrays. 

Ultra-high-resolution radio observations provide a route to probing the inner- 
most regions of ULIRGs (Smith et al., 1998; Carilli and Taylor, 2000). By 
detected the diffuse emission and multiple point-like radio sources, expected 
from multiple supernova remnants, rather than a single point-like core and 
accompanying jet structures expected from an AGN, these observations sug- 
gest that high-mass star formation contributes at least a significant part of 
the luminosity of the ULIRGs Arp220 and Mrk273. 

It is interesting to note that the observed correlation between the inferred mass 
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of the black holes in the centers of galaxies and the stellar velocity dispersion 
of the surrounding galactic bulges, in which most of the stars in the Universe 
reside (Fukugita et al., 1999), might inform this discussion (Magorrian et al., 
1998; Ferrarese and Merritt, 2000; Gebhardt et al, 2000). The mass of the 
bulge appears to exceed that of the black hole by a factor of about 200. When 
hydrogen is processed in stellar nucleosynthesis, the mass-energy conversion 
efficiency is about 0.007e*, where e*(~ 0.4) is the fraction of hydrogen burned 
in high-mass stars. When mass is accreted onto a black hole, the mass-energy 
conversion efficiency is expected to be about 0.1€bh, with 6bh ~ 1 with the 
definition above. If accretion and nucleosynthesis were to generate the same 
amount of energy during the formation of a galaxy, then the ratio of mass 
contained in both processed stars and stellar remnants to that of a super- 
massive black hole is expected to be about O.leBH/0.0076*. For e* = 0.4 and 
e B H = 1, this ratio is about 36. As a mass ratio of about 200 is observed, this 
implies that a greater amount of energy, by a factor of about 6, is generated 
by high-mass star-formation activity than by gravitational accretion. 

If the bulge-to-black-hole mass ratio is in fact greater than 200, then either the 
factor by which star formation dominates will exceed 6, or the accretion must 
have been more than 10% efficient; that is e B H > 1- If low-efficiency accretion 
dominates the process of the build up of mass in the central black hole, then 
less than 1 part in 7 of the luminosity generated during galaxy formation will 
be attributed to accretion as compared with high-mass star formation. 

A greater amount of energy generated by star formation as compared with 
accretion processes appears to be favored by these circumstantial arguments. 



2.2 Continuum emission from dust 

The dust emission process is thermal, with dust grains emitting a modified 
blackbody spectrum. Grains of interstellar dust, distributed throughout the 
ISM of a galaxy, are heated to temperatures between about 20 and 200 K, 
depending on the spectrum and intensity of the interstellar radiation field 
(ISRF), and the size and optical properties of the grains. Higher dust tem- 
peratures can be produced close to a powerful source of radiation, with dust 
subliming at temperatures of order 2000 K. Very small grains can be heated 
far above their equilibrium temperatures by absorbing hard-UV photons (see 
Draine and Li, 2001). Lower dust temperatures, always exceeding the CMB 
temperature, are possible in opaque regions of the ISM that are shielded from 
intense heating, in the intergalactic medium or in regions with an intrinsically 
weak ISRF. Unless dust is heated by the ISRF in addition to the CMB the 
galaxy will not be detectable. We now consider the properties of the dust 
emission that are relevant to observations of high-redshift galaxies. 
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2.2.1 The emission spectrum, dust mass and temperature 

The minimum parameters necessary to describe the emission from dust grains 
are a temperature Td and a form of the emissivity function e v . In any galaxy 
there will be a distribution of dust temperatures, reflecting the different na- 
ture and environment of each grain. It is useful to use Td to describe the 
coolest grains that contribute significantly to the energy output of a galaxy 
when discussing submm observations. In most cases, spatially and spectrally 
resolved images of galaxies are not available, and so it is reasonable to as- 
sume a volume-averaged description of the emissivity function as a function of 
frequency u, e„ oc iA Values of (3 in the range 1-2 are usually assumed. Scat- 
tering theory predicts that (3 — > 2 at low frequencies, while a value f3 ~ 1 at 
high frequencies matches the general trend of the interstellar extinction curve 
that describes the properties of absorption of optical and UV radiation by 
the ISM (see Calzetti et al., 2000 and Section 2 of the review by Franceschini, 
2002). 

The simplest form of the emission spectrum/SED, f v is given by assuming that 
f v oc e u B u , in which B v is the Planck function (2kT d u 2 /c 2 in the Rayleigh- 
Jeans limit, in units of Wm~ 2 Hz -1 sr _1 ). This assumes that the emitting 
source is optically thin. For fitting spectra of galaxies found in deep submm 
surveys, we assume the simple t v B v function to describe the SED. Dunne et 
al. (2000) and Dunne and Eales (2001) also use this functional form to fit the 
observed submm spectra of low-redshift galaxies. At the expense of adding 
another parameter to describe the SED, there is some physical motivation for 
a SED that includes an optical depth term, 

/„ oc [1 - exp(-T„)]B v , (1) 

where r v is the frequency-dependent optical depth of the cloud, and is a multi- 
ple of e v . This equation tends to the simpler e v B v function at long wavelengths, 
and is assumed by, for example, Benford et al. (1999), Omont et al. (2001), 
Priddey and McMahon (2001) and Isaak et al. (2002), whose submm data 
for high-redshift AGNs tends to correspond to rest-frame frequencies that are 
relatively close to the peak of the SED. The extra parameter required to re- 
late t v and e v can be defined as the frequency at which t v — 1 and the cloud 
becomes optically thick. If the opacity near a wavelength of 100 /im is large, 
then the form of the peak of the SED tends to that of a blackbody spectrum. 
This suppresses the emission near to the SED peak relative to the emission in 
the Rayleigh- Jeans regime, and so this functional form provides a good fit to 
a set of submm and far-IR data with a higher value of Td as compared with 
the e v B v function, usually by about 10-20%. However, because most observed 
SEDs for high-redshift galaxies have fewer than four data points (see Fig. 2), 
the difference is unlikely to be very significant. 
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It is reasonable to assume that the mid-IR SED can be smoothly interpolated 
from a modified blackbody function at low frequencies to a power-law f v oc u a 
in the mid-IR waveband on the high-frequency side of the spectral peak, in 
order to prevent the high-frequency SED from falling exponentially with a 
Wien spectrum. Hotter components of dust, emitting at shorter wavelengths, 
and ultimately stellar emission in the near-IR waveband, are certain to be 
present to reduce the steepness of the SED in the Wien regime. That an 
exponential Wien spectrum is inappropriate can be seen from the well-defined 
power-law mid-IR SEDs of Arp220 and Mrk231 shown in Fig. 2. 

It is not always necessary to relate the SED f v and luminosity L of a galaxy 
to the mass of dust that it contains; this can of worms can remain closed 
by normalizing /„ in a self-consistent way. However, if a dust mass is required, 
perhaps in order to estimate the metal content of the ISM, and so provide 
information about the integrated star- format ion activity in the galaxy at ear- 
lier times (Hughes et al., 1997; Omont et al., 2001), then it is conventional to 
define a frequency-dependent mass-absorption coefficient k u (Draine and Lee, 
1984; with units of m 2 kg _1 ), which is proportional to e u . k v is the 'effective 
area' for blackbody emission by a certain mass of dust, 



Values of k v at a conventional frequency of around 1 mm are in the range 0.04- 
0.15 m 2 kg -1 (Hughes, 1996). Recent comparisons of optical extinction and 
submm emission from partially resolved edge-on spiral galaxies have tended 
to give values of 0.05-0.4 m 2 kg -1 (see Fig. 4 of Alton et al., 2001). Domingue 
et al. (1999) derive 0.09m 2 kg _1 from similar far-IR, optical and submm data. 
Dunne et al. (2000) adopt a value of 0.077 m 2 kg -1 . Note that there is at least 
a factor of 3 uncertainty in these conversion factors. 

An alternative dimensionless function Q u (Hildebrand, 1983) is sometimes 
used, which includes information about the mass/volume and surface area of 
a typical grain. If grains are assumed to be spherical (a big if), with bulk 
density p, radius a, and an emissive cross section ira 2 , then k v = 3Q u /4ap. 
Q U B U is the effective emissivity function describing the energy flux from unit 
area of the dust grain surface. However, dust grains are more likely to be 
irregular in shape, possibly colloidal or in the form of whiskers. In that case, 
the emissivity per unit mass would be increased, and the dust mass associated 
with a fixed luminosity would be overestimated. 

This geometrical uncertainty will inevitably result in uncertainty about the 
mass of dust. Hence, dust masses quoted in papers must be treated with cau- 
tion, and may be best used as a comparative measure to distinguish galaxies. 
In general, we will avoid quoting dust masses, as this is unlikely to provide a 
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reliable physical measure of the properties of galaxies until detailed resolved 
images are available, which is likely to require observations with the ALMA 
interferometer. This will be a recurring theme: observations with excellent sen- 
sitivity and spatial resolution using a large interferometer will resolve many 
of the questions raised throughout the paper. 

Working from submm data, it is also difficult to assess the dust mass of a 
galaxy, even subject to the caveats above, without knowing its dust temper- 
ature. In the Rayleigh- Jeans spectral regime, the flux density from a galaxy 
S v oc v 2+/3 M d T d . If T d is uncertain to within a factor, then M d is uncertain 
to within the same factor. The dust mass is at least easier to estimate from a 
single long- wavelength observation than the luminosity L. As L oc M d T d +/3 , or 
equivalently L oc S u T^ +f3 , an uncertainty in T d corresponds to a proportionally 
much larger uncertainty in the inferred value of L. 

However, even if the dust mass can be determined reliably at low redshifts, it 
remains unclear whether the same procedure can be applied to determine the 
dust mass in more luminous and more distant systems. In order to determine 
the dust properties of high-redshift galaxies, data of the same quality that 
has been obtained for nearby galaxies is required. High-frequency submm/far- 
IR observations are necessary to provide information about the rest-frame 
frequency of the peak of the SED for a high-redshift galaxy. 

Given the current lack of resolved images of distant galaxies in the submm and 
far-IR wavebands, it is important to neither over-parametrize the descriptions 
nor overinterpret the results of observations of their SEDs. When spatially- 
resolved, high-spectral resolution images are available, building on existing 
interferometric images of low-redshift dusty galaxies (Downes and Solomon, 
1998; Sakamoto et al., 1999; Wilson et al., 2000), it should be possible to study 
the radiative transfer from sites of intense star formation and AGN in these 
geometrically complex opaque galaxies (see Ivison et al., 2000a, 2001). Models 
of the SEDs of dust-enshrouded AGN at different viewing angles have been 
developed by Granato et al. (1996), while star-forming regions embedded in a 
disk geometry have been analyzed by Devriendt et al. (1999). More powerful 
and efficient radiative-transfer codes are being developed (for example Abel 
et al., 1999), and it should be practical to develop detailed models of the 
appearance of galaxies with realistic geometries to account for future, high- 
resolution multi-band submm images. 

At present, we prefer to use a few simple parameters — a, f3 and T d — to describe 
the essential features of the SEDs of dusty galaxies. Although such a model 
can encapsulate only a small part of the true complexity of the astrophysics 
in a galaxy, it can account for the existing SED data for a wide variety of 
dusty galaxies. A simple parametrization is preferable to a more baroque, 
and necessarily at present unconstrained, combination of geometry, dust mass 
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and temperature. In the following section we list plausible values of our SED 
parameters and discuss the associated degeneracies in fitted values. 



2.3 The observed SEDs of dusty galaxies 

Information about the submm SEDs of galaxies has been gathered from tar- 
geted mm and submm observations of samples of low-redshift far-IR-selected 
galaxies from the IRAS catalog, (Andreani and Franceschini, 1996; Dunne et 
al., 2000; Lisenfeld et al., 2000; Dunne and Eales, 2001), and from far-IR and 
submm observations of high-redshift galaxies (see Fig. 2). The most extensive 
local survey (SLUGS; Dunne et al., 2000) consists of 850-/zm SCUBA obser- 
vations of 104 galaxies selected from the low-redshift IRAS Bright Galaxy 
Sample (BGS; Soifer et al., 1987). After fitting single-temperature e u B v SEDs 
to the galaxies, Dunne et al. found that (3 — 1.3 ± 0.2 and T d = 38 ± 3K 
described the sample as a whole, with a natural dispersion in the properties 
from galaxy to galaxy. This 7iL45'-selected sample could be biased against less 
dusty galaxies. Dunne et al. are currently addressing this issue by observing a 
complementary sample of £>-band selected low-redshift galaxies, which should 
be representative of optically luminous low-redshift galaxies as a whole. 

Note, however, that when fitting only a few datapoints, there is a significant 
correlation between values of (5 and Td that can account for the data (left panel 
of Fig. 3). This can lead to ambiguity in the results, further emphasizing the 
difficulty in associating the dust mass or temperature inferred from a galaxy 
SED with the real physical properties of the galaxy. 

The addition of 450-yum data for 19 of the 104 galaxies in the SLUGS sample 
(Dunne and Eales, 2001), tends to split the galaxy SEDs into two categories: 
those that retain a definite 40-K spectrum after including the 450-/im data, 
and those for which cooler single-temperature SEDs, more similar to the SEDs 
of normal spiral galaxies, then provide a better fit. The first group are typically 
the more luminous galaxies in the sample, while the second includes 3 of the 5 
lowest luminosity galaxies from the sample. Dunne and Eales (2001) propose 
a two-temperature model to account for the changes in light of the new 450- 
[xm data; however, a cooler single-temperature model with a larger value of f3 
provides a fit of similar quality. The results for one of the most significantly 
different fits is shown in the right-hand panel of Fig. 3. With the addition 
of the 450-/iin data, the nature of the SEDs of low-redshift, low-luminosity 
galaxies become more diverse. However, the more luminous galaxies, which 
are likely to be the most similar to typical high-redshift submm galaxies, are 
still described reasonably well by the original Dunne et al. (2000) 38-K SED. 

An alternative approach is to determine an SED that can describe the ob- 
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Fig. 3. An illustration of some of the issues involved in describing the SEDs of dusty 
galaxies. On the left is a probability contour plot that shows the 0.5, 5 x 1CT 3 and 
5 x 10~ 5 probability contours for a fit to an SED model defined by the variable 
parameters (3 and Td with a fixed value of a = —1.95, taking into account four 
SED datapoints for the galaxy NGC 958 as shown in the right-hand panel (Dunne 
and Eales, 2001). Note that Uy = 10~ 26 Wm~ 2 Hz" 1 . Note that there is a very 
significant degeneracy in the fitted parameters. Adding additional data points with 
small errors close to the peak of the SED at 200 /im reduces the extent of the 
probability contours by about 50%, but they remain elongated in the same direction. 
Note that (3 > 2 is not expected physically. On the right the data are compared 
with fitted single-temperature SEDs. The solid line is the best fit to the data. 
The dashed lines correspond to SEDs from the ends of the probability 'banana' 
shown in the left-hand panel. Note that without the 450-^/m point, the thick dashed 
curve describes the best-fit SED, which is defined by a significantly greater dust 
temperature. This SED is similar to that of a typical luminous IR galaxy, whereas 
the best fitting model with all four data points is much more like the SED of the 
Milky Way. Note that the shift in the best-fit model on adding 450-/im data is 
generally less significant than in this case. 

served flux density distribution of galaxies in the far-IR and submm wave- 
bands, which are sensitive to galaxies at low, moderate and high redshifts 
(Blain et al., 1999b; Trentham et al., 1999; Barnard and Blain, 2002). Using 
the e v B v functional form, values of (5 ~ 1.5 and Ti ~ 40 K are required to 
provide a good description of the data, rather similar to the values derived 
for temperatures of individual low-redshift luminous dusty galaxies in Dunne 
et al. (2000) and Lisenfeld et al. (2000), and for both the small number of 
high-redshift submm-selected galaxies with known redshifts and mid-IR spec- 
tral constraints (Ivison et al., 1998a, 2000a) and typical high-redshift QSOs 
(for example Benford et al., 1999). These temperatures are significantly less 
than those determined for the most extreme high-redshift galaxies (Lewis et 
al., 1998), and significantly greater than the Td = 17 K inferred from the 
maps of the Milky Way made using the all-sky survey from the FIRAS in- 
strument on the Cosmic Background Explorer ( COBE) satellite in the early 
1990's (Reach et al., 1995). Note that there are examples of moderate-redshift 
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infrared-selected galaxies with both hotter and colder typical dust tempera- 
tures than 40 K: see Deane and Trentham (2001) and Chapman et al. (2002d) 
respectively. At present it seems likely that a 40-K dust temperature is a 
reasonable assumption for high-redshift submm-selected galaxies. 

Inevitably, however, there will be a population of hotter high-redshift galaxies 
(Wilman et al., 2000; Trentham and Blain, 2001). These galaxies would be 
underrepresented in existing submm surveys, but may make a significant con- 
tribution to the 240-/im background radiation intensity (Blain and Phillips, 
2002). Further observational information to test the assumption of a 40-K 
dust temperature is keenly awaited. As we discuss below, in Section 2.6, the 
assumed dust temperature has a significant effect on the selection function of 
submm galaxy surveys, and on the properties that are inferred for the galaxies 
that are found in these surveys. 

2.4 Line emission 

Emission from molecular rotation and atomic fine-structure transition lines 
can be used to diagnose physical conditions within molecular clouds and pho- 
todissociation regions, and to trace out the velocity structure within. Some 
lines, such as those from CS, HCN and HCO + are excited only in high-density 
gas, while others, including the most abundant polar species CO, trace more 
typical regions in the ISM. 

Studies of many emission lines from molecular cloud regions in nearby galaxies 
are possible using existing mm and submm- wave telescopes (Wilson et al., 
2000; Heifer, 2000). However, for more distant galaxies only CO lines have 
so far been detected in significant numbers, almost exclusively from galaxies 
which have been subject to strong gravitational lensing by foreground galaxies 
(see the summary in Combes et al., 1999). These observations are useful for 
deriving physical conditions within the sources, especially if multiple lines 
are detected (as in the case of APM 08279+5255; Downes et al. 1999b). The 
improved capabilities of the forthcoming mm/submm interferometer arrays — 
SMA, upgrades to the IRAM Plateau du Bure interferometer (PdBI), and the 
Combined Array for Research in Millimeter-wave Astronomy (CARMA)Q- 
and and ultimately the dramatically increased sensitivity of ALMA, will make 
high-redshift lines much easier to observe over the next decade (Combes et al., 
1999; Blain et al., 2000b). 

One of the most important uses of CO-line observations of distant submm 
galaxies found in continuum surveys is their ability to confirm an identification 
absolutely, by tying together an optical and submm redshift at the position 

7 http://www.mmarray.org/ 
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of the galaxy. So far this has been achieved for only three submm galaxies 
(Frayer et al, 1998, 1999; Kneib et al., 2002: see Figs. 14, 18 and 19). In 
principle, these observations could be made for all continuum-selected galaxies. 
The difficulty is the narrow fractional bandwidth available for receivers and 
correlators. Even at the relatively low frequency of 90 GHz, the redshift of 
the target must be known to better than 0.5% to ensure that a 300 km s -1 
wide CO line, typical of a massive galaxy, with a width equivalent to 0.1% 
in redshift falls entirely within a 1-GHz band. Future cm-, mm-, and submm- 
wave instruments with wider bandwidths will significantly assist the the search 
for redshifts using molecular lines. Specially designed low-resolution, ultra- 
wideband dispersive spectrometers covering many tens of GHz simultaneously 
on single-antenna mm-wave telescopes also promise to provide redshifts for 
submm galaxies (Glenn, 2001). 

A complementary search for redshifted cm-wave OH megamaser emission to 
pinpoint the redshifts and positions of ultraluminous high-redshift galaxies 
could be possible using radio telescopes (Townsend et al., 2001). However, 
there are very stringent requirements on the acceptable level of radio fre- 
quency interference from terrestrial and satellite communications. Observa- 
tions of low-redshift megamasers are described by Darling and Giovanelli 
(2001). Megamaser emission at high redshifts is discussed by Briggs (1999) in 
the context of the proposed Square Kilometer Array (SKA) meter /centimeter- 
wave radio telescope. If it can operate at frequencies of several 10's of GHz, 
then the SKA is also likely to be an efficient detector of low-excitation high- 
redshift CO lines (Carilli and Blain, 2002). 



2.4-1 Line emission contribution to continuum detections 

An interesting feature of the CO line emission from low-redshift galaxies is 
that lines can lie in the passbands of continuum instruments, and could con- 
tribute to the continuum flux inferred. For low-redshift galaxies, the 345- GHz 
CO (3 — > 2) line lies within the 850-/xm atmospheric window, while the 691- 
GHz CO(6 -> 5) and 230-GHz CO(2 -> 1) lines lie in the 450-//m and 1.25-mm 
windows, respectively. 

Assuming a reasonable template spectrum (Blain et al., 2000b), the equivalent 
width in frequency of the CO(3 — > 2) line is 7.4 GHz. The passband of the 
current SCUBA 850-/im (353-GHz) filter is about 120 /im (50 GHz) wide, and 
so about 15% of the measured continuum flux density of a low-redshift galaxy 
in the 850- /im channel is likely to be from the CO line. The high-frequency 
SCUBA passband in the 450-/im atmospheric window is 75 GHz wide, while 
the equivalent width of the CO(6 — > 5) transition is 3.3 GHz. Hence, a smaller 
5% contribution to the continuum flux density from the line is expected at 
450 fxm. The CO(2 — > 1) line has an expected equivalent width of 9.2 GHz, 
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while the wide MAMBO passband has half-power points at 210 and 290 GHz. 
Contamination of the flux densities detected by MAMBO by about 10% may 
thus be expected. 

The largest of these correction factors is comparable to the calibration un- 
certainty in submm-wave observations, and could be relevant to the detailed 
interpretation of low-redshift observations. For example, the presence of the 
CO (3 — > 2) line in the 850-/zm window would shift the inferred continuum 
emissivity spectral index (3 in the SLUGS survey from 1.3 to 1.52. At high 
redshifts, any corrections are likely to be less significant, both because the 
relatively bright CO (3 — > 2) line redshifts out of the 850-yum passband, and 
the equivalent width of lines in frequency space decreases as (1 + z)^ 1 . 

Although the contribution to measured submm-wave flux densities from line 
emission could be significant at the level of order 10%, only a small fraction of 
the bolometric luminosity from galaxies is detected in the submm waveband. 
More than 99% of the bolometric luminosity still appears in the continuum, 
predominantly at shorter far-IR wavelengths. 



2.5 The observability of high-redshift dusty galaxies 



The detectable flux density at an observed frequency v from a galaxy with 
bolometric luminosity L at redshift z with an intrinsic SED f u , 

C _ 1 + Z T /"(!+*) (o\ 



where D L is the luminosity distance to redshift z (for example Peebles, 1993). 

The key feature that makes submm-wave observations of distant galaxies in- 
teresting is the ability to sample the SED of a target galaxy at wavelengths for 
which the SED is a strongly increasing function of frequency (Fig. 2). This en- 
sures that distant galaxies are observed at a rest-frame wavelength closer to the 
peak of their SED. There is thus a strong, negative K correction, which leads 
to high-redshift galaxies being relatively easy to detect at submm wavelengths 
as compared with their low-redshift counterparts. This effect is illustrated in 
Fig. 4 for the template SED from Blain et al. (1999b) shown in Fig. 2. The 
strong K-correction effect applies at wavelengths longer than about 250 fxm. 
At these wavelengths the flux density from galaxies at z > 1 ceases to decline 
with the inverse square of distance, but instead remains approximately con- 
stant with increasing redshift. A window is thus opened to the detection of all 
galaxies with similar SEDs at redshifts up to z ~ 10-20. The effect is more 
pronounced at longer wavelengths: in the mm waveband more distant galax- 
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Fig. 4. The predicted flux density of a dusty galaxy as a function of redshift in various 
submm atmospheric windows, and at shorter wavelengths that will be probed by 
forthcoming space missions. Note the powerful K correction in the mm and submm 
wavebands at wavelengths longer than about 250 fim, which yields a flux density 
that is almost independent of redshift. The template spectrum is chosen to reproduce 
the typical properties of distant submm-selected galaxies (Fig. 2). Subtle effects due 
to the additional heating of dust by the CMB, and fine details of the radio SED of 
galaxies are not included; these effects are illustrated in Fig. 8. 

ies are expected to produce greater flux densities than their more proximate 
counterparts. 

Note that both the radio and optical flux-density-redshift relations decline 
steeply with increasing redshift, and so high- redshift galaxies are not selected 
preferentially in those wavebands. The advantage that faint radio and optical 
galaxy surveys have over submm surveys comes from the complementary probe 
of astrophysical signatures, and the combination of greater fields of view and 
finer angular resolution. 

A submm telescope that is sufficiently sensitive to detect a certain class of 
galaxy at redshift z ~ 0.5, can detect any similar galaxies out to a redshift 
z ~ 10 (Blain and Longair, 1993a). Note, however, that surveys to exploit this 
unusual K correction are not immune to selection effects. The K correction 



19 



can also only be exploited at redshifts for which sufficient heavy elements are 
present in the ISM of the target galaxy to form enough dust to reprocess opti- 
cal radiation. Nor does the K correction effect overcome cosmological surface 
brightness dimming for progressively more distant submm galaxies: the nor- 
mal (1 + z)~ A reduction in surface brightness still applies; however, it is not 
expected to become significant until redshifts in excess of about 5. Because 
submm-wave telescopes do not yet resolve distant galaxies, this effect cannot 
be observed at present. It may provide an opportunity to estimate redshifts 
for the most distant submm-selected galaxies when they can be resolved using 
ALMA. 



2.6 Submm-wave selection effects 

Deep submm-wave observations image the high-redshift Universe with very 
little contamination from low-redshift galaxies, and can potentially find a pop- 
ulation of galaxies that is quite different to those detected in conventional deep 
optical surveys, and which could be undetectable in these surveys. The comple- 
mentarity of submm and optical observations is illustrated by the very limited 
overlap between galaxies detected in the deep submm-optical image shown in 
Fig. 1. However, submm surveys are certainly subject to selection effects. In 
Fig. 5 the flux-density-redshift relation for a submm-luminous galaxy with a 
fixed bolometric luminosity is presented as a function of its SED parameters 
- T d , a and (3. The relatively minor effects of different assumed cosmological 
models are also shown. Changing the dust temperature has the greatest ef- 
fect. The inferred luminosity of a dusty galaxy for a fixed observed submm flux 
density goes up by a factor of 10 if the dust temperature is doubled, at all but 
the very highest redshifts. There is thus a significant potential bias in submm 
surveys against the detection of galaxies with hotter dust temperatures for a 
given bolometric luminosity. 

This effect was noted by Eales et al. (1999), when investigating the evolution of 
galaxies in the context of the results of deep SCUBA surveys. They suggested 
that the submm galaxies may be cooler than the temperatures of about 60 K 
usually assumed, and so their significance as a population of strongly evolving 
high-redshift galaxies may have been overestimated. 

As discussed in Section 2.3, a cooler dust temperature of 40 K is compatible 
with observations of the SEDs of individual submm galaxies with confirmed 
redshifts detected in submm surveys (Ivison et al., 1998a, 2000a) and with the 
results of targeted observations of luminous low-redshift IRAS galaxies and 
high-redshift QSOs. If this temperature is assumed, then the inferences about 
galaxy evolution made from the results of submm surveys (Blain et al., 1999b, 
c; Eales et al., 2000; Smail et al., 2002) should be reliable. However, until a 
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Fig. 5. Flux-density-redshift relations illustrating some of the uncertainties that 
apply to the interpretation of submm and far-IR surveys. In the top-left panel the 
relatively small effects of changing the world model parameters are shown. Note that 
changes to the volume element and flux density received counteract each other, and 
so the chosen cosmology has a small effect on the interpretation of surveys. The most 
significant effect — that of changing the dust temperature — is shown in the top-right 
panel. At 175 /um for galaxies at moderate redshifts, the effect of temperature is 
rather small. However, at 850 ^m, the effects are very significant, and must be 
remembered when interpreting the results of 850-^m observations: doubling the 
dust temperature corresponds to increasing the luminosity associated with a given 
flux density by a factor of about 10. The less significant effects of changing the dust 
emissivity index (3 or the mid-IR spectral index a are shown in the bottom-left and 
bottom-right panels respectively. 

large sample of submm galaxies with redshifts and multi-waveband SEDs is 
available, the possibility that a cold or hot population of high-redshift dusty 
galaxies could be missing from or misidentified in submm surveys cannot be 
ruled out (Eales et al., 1999; Blain and Phillips, 2002). The possible effects 
on inferred luminosities of different forms of the SED shown in Fig. 5 need 
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to be taken seriously, especially when describing the properties of individual 
galaxies selected in submm surveys. 

There is little reliable evidence for a systematic relationship between dust tem- 
perature and redshift. Observations of low-redshift IRAS galaxies (Andreani 
and Franceschini, 1996; Dunne et al., 2000), indicate that any variation of 
dust temperature with luminosity appears to be gradual. However, there is 
evidence for a significant and systematic change in the temperature of dusty 
galaxies with a wider range of luminosities, from about 20 K for low-redshift 
spirals (Reach et al., 1995; Alton et al., 2000; Dunne and Eales, 2001) to about 
40 K for more luminous objects typical of the galaxies detected in the IRAS 
survey. Temperatures of up to 110 K are found for some extremely luminous 
high-redshift galaxies (Lewis et al., 1998). 

We stress that there could be a significant selection effect in submm surveys 
that depends on the range of dust temperatures in the source population. 
The importance of such an effect can be quantified once a complete redshift 
distribution is available for a submm-selected galaxy sample. 



2.7 Deep submm-wave surveys 

Images of the redshift z = 0.25 cluster of galaxies Abelll835 in both the op- 
tical and submm wavebands were compared in Fig. 1. This image provides a 
realistic impression of the appearance of deep optical and submm images of 
the sky. Note that the only relatively low-redshift or cluster member galaxy 
that contributes any submm-wave flux is the central cD galaxy: the other clus- 
ter galaxies are either quiescent, neither forming stars nor heating dust, or are 
insufficiently luminous to be detectable at 850 /zm using SCUBA. Background 
galaxies at much greater redshift, which have faint optical counterparts as 
compared with the cluster member galaxies, dominate the image. This is a 
direct visual demonstration of the strong bias towards the detection of distant 
galaxies in submm-wave surveys that was illustrated in Fig. 4. These back- 
ground galaxies are magnified by a factor of order 2-3 due to the gravitational 
lensing potential of the foreground cluster over the full extent of the image. 
The effects of gravitational lensing can be determined using accurate mod- 
els of the cluster potential, that are constrained with the help of data from 
HST images and spectroscopic redshifts for multiply-imaged optically-selected 
galaxies. The uncertainty in the results is comparable to the uncertainty in 
the calibration of the submm images. 

As shown in Fig. 1, existing deep submm images are much less visually stimu- 
lating than deep optical images, because their angular resolution is not suffi- 
cient to image the internal structure in distant galaxies. The limited resolution 
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also imposes a confusion limit to the depth for submm surveys, at which the 
noise level is dominated not by atmospheric or instrumental noise but by the 
telescope resolution blurring together signals from faint unresolved galaxies. It 
takes about 50 h of integration using SCUBA to reach the practical confusion 
limit in a single field. Confusion is discussed in more detail in Section 3.1. 

There are a variety of published results from deep submm galaxy surveys. Sur- 
veys aim to detect high-redshift galaxies exploiting the powerful K-correction 
effect in the submm waveband. Over 500 arcmin 2 of blank sky has been sur- 
veyed using SCUBA by several groups (Barger et al., 1998; Hughes et al., 1998; 
Barger et al, 1999a; Eales et al, 1999, 2000; Borys et al, 2002; Fox et al, 2002; 
Scott et al., 2002; Webb et al., 2002a). These range from an extremely deep 
survey in the area of the Hubble Deep Field-North (HDF-N; Williams et al. 
1995) by Hughes et al. (1998) searching for the faintest detectable populations 
of submm galaxies, to wider-field shallower surveys to detect brighter sources 
that might be easier to follow-up and could be used to trace large-scale struc- 
ture (Borys et al. 2002; Scott et al. 2002). About 30 5-arcmin 2 lensed cluster 
fields have been imaged using SCUBA (Smail et al., 1997, 2002; Chapman et 
al., 2002a; Cowie et al., 2002; Kraiberg Knudsen et al., 2001; van der Werf and 
Kraiberg Knudsen, 2001), to various RMS depths between 0.5 and 8mJy.[^] 
By exploiting the magnification effect of gravitational lensing, which extends 
over fields several arcminutes across, due to rich clusters of galaxies at mod- 
erate redshifts, the population of distant galaxies in the source plane behind 
the lensing cluster can be probed to greater depths than is possible in a blank 
field (Blain, 1998). 

The detection rate of galaxies using SCUBA based on published papers ap- 
pears to have declined over time since 1998. In significant part this is due to 
the absence of the sustained excellent observing conditions on Mauna Kea that 
were experienced during the El Nino winter of 1997-1998, just after SCUBA 
was commissioned. 

A wide variety of complementary, and sometimes overlapping surveys have 
been made using MAMBO at the IRAM 30-m telescope (Bertoldi et al, 2000, 
2001; Carilli et al., 2001), during several winters. These fields include the 
cluster Abell2125 and the ESO-NTT Deep Field (Arnouts et al., 1999). 

A compilation of the results from all the SCUBA and MAMBO surveys is 
presented in Figs. 9 and 10. A full summary of deep projects that have been 
undertaken or are underway can be found in Ivison (2001). In addition, larger 
shallower submm surveys of the Galaxy (Pierce-Price et al., 2000), and per- 

8 Uy = lO^Wm^Hz- 1 

9 A JCMT project (Vicki Barnard et al.) is currently searching for candidate high- 
redshift galaxies detected in wide-field SCUBA images of star-forming regions in 
the Milky Way (Barnard et al., 2002). 
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haps the CMB images obtained using the BOOMERANG balloon-borne ex- 
periment (Masi et al., 2001), can be used to search for brighter submm-wave 
galaxies. 

2.8 Submm observations of known high-redshift galaxies and QSOs 

The advent of SCUBA and MAMBO has also provided the opportunity to 
study the submm properties of large samples of interesting high-redshift galax- 
ies, including almost all types of previously known distant galaxies. Isolated 
detections of high-redshift AGN-powered radio galaxies and QSOs were made 
in the mid 1990 's using single-element bolometer detectors (for example Dun- 
lop et al., 1994; Isaak et al., 1994); however, the compilation of statistical sam- 
ples, and the secure rejection of contamination from fluctuating atmospheric 
noise have only been possible more recently, using SCUBA and MAMBO, and 
the 350-/im one-dimensional bolometer array SHARC at the 10.4-m aperture 
Caltech Submillimeter Observatory (CSO) on Mauna Kea. A key advantage 
of observing these sources is that both their redshifts and some of their astro- 
physical properties are already known, in contrast with the submm-selected 
galaxies discovered in blank-field surveys. Some of the targeted galaxies — 
very faint non-AGN radio galaxies, mid-IR-selected ISO galaxies, and X-ray 
selected AGNs — have only been detected very recently. As the relationship 
between these populations of galaxies and submm-selected galaxies is still un- 
clear, many of the limits will be discussed in the context of following up submm 
surveys in Section 4. 

Targeted surveys include a search for submm-wave continuum emission from 
high-redshift AGN-powered radio galaxies (Archibald et al., 2001), and obser- 
vations of various samples of optically-selected QSOs (for example Carilli et 
al., 2001; Isaak et al., 2002). In these observations a single bolometer is aimed 
at the position of the target. While this does not lead to a fully-sampled im- 
age of the sky, it provides a more rapid measurement of the flux density at a 
chosen position. The results have been the detection and characterization of 
the dust emission spectra for a range of luminous high-redshift galaxies and 
QSOs, including APM 08279+5255 (Lewis et al, 1998), the galaxy with the 
greatest apparent luminosity in the Universe. Barvainis and Ivison (2002) have 
targeted all the known galaxies magnified into multiple images by the gravita- 
tional lensing effect of foreground galaxies from the CASTLES gravitational 
lens imaging project,^] significantly expanding the list of high-redshift galax- 
ies magnified by a foreground mass concentration with a submm detection. 
The SEDs of several of these galaxies are shown in Fig. 2. 

Archibald et al. (2001) find evidence for significant evolution in the properties 
10 http:/ /cfa- www. harvard.edu/castles. 



24 



of dust emission with increasing redshift in a carefully selected sample of AGN 
radio galaxies, whose radio properties were chosen to be almost independent 
of the redshift of the observed galaxy. The results perhaps indicate that more 
intense star- format ion activity, as traced by the submm emission, takes place 
alongside the radio source activity at higher redshifts, and so provide a possible 
clue to the formation and evolution of the massive elliptical galaxies thought 
to host radio galaxies. Hughes et al. (1997), Ivison et al. (1998b) and Omont 
et al. (2001) discuss the consequences of finding large masses of dust at high 
redshifts, in terms of the limited cosmic time available for the formation of the 
stars required to produce the metals and dust required to generate sufficiently 
intense submm emission from the host galaxy. 

The radio galaxies detected by Archibald et al. (2001) in pointed single- 
bolometer SCUBA observations were followed up by imaging observations of 
the surrounding 5-arcmin 2 fields, to search for submm-loud companions. Ivi- 
son et al. (2000b) found that the surface density of submm galaxies in some 
of these fields is about an order of magnitude greater than that in a typical 
blank field, indicating a significant overdensity of sources. This is likely due to 
some radio galaxies being found in high-density regions of biased high-redshift 
galaxy formation, which are possibly 'protoclusters' — rich clusters of galaxies 
in the process of formation. 

A similar targeted approach has been taken to try to detect submm-wave 
emission from optically-selected LBGs at redshifts between 2.5 and 4.5 (Steidel 
et al., 1999). The Lyman-break technique (Steidel et al., 1996) detects the 
restframe 91.2-nm neutral hydrogen absorption break in the SED of a galaxy as 
it passes through several broad-band filters. Large samples of candidate LBGs 
can be gathered using multi-color optical images from 4-m class telescopes. 
The efficiency of the selection method is of order 70% after spectroscopic 
confirmation of the candidates using 8/10-m class telescopes. The LBGs are 
the largest sample of spectroscopically confirmed high-redshift galaxies, with 
a well-defined luminosity function (Adelberger and Steidel, 2000), a surface 
density of order lOarcmin" 2 , and inferred star- formation rates between 1 and 
10M Q yr _1 . They appear to be typical of the population of distant galaxies, 
and their spectra provide useful astrophysical information. 

Observing the LBGs at submm wavelengths is an important goal, clS £111 duC- 
curate determination of their submm-wave properties will investigate the link 
(if any) between the large well-studied LBG sample and the more enigmatic 
submm galaxy population (Blain et al., 1999c; Lilly et al., 1999; Adelberger 
and Steidel, 2000; Granato et al., 2001). At present, the typical very faint 
limits to the optical counterparts of the subset of submm galaxies with accu- 
rate positions (Smail et al., 1998a, 2002; Downes et al., 1999; Dannerbauer et 
al., 2002), and the detection of Extremely Red Object (ERO) galaxies (with 
R — K > 6) as counterparts to a significant fraction of submm galaxies (Smail 
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et al, 1999, 2002; Gear et al., 2000; Frayer et al, 2002; Ivison et al, 2001; 
Lutz et al., 2001), argue against a large overlap between the two populations. 
The direct submm detection of LBGs using SCUBA has been largely unsuc- 
cessful at the 0.5-mJy RMS level: a single galaxy out of 16 was detected by 
Chapman et al. (2000; 2002c), while Webb et al. (2002b) describe a low signif- 
icance of overlap between LBGs and SCUBA galaxies in a wide-field survey. 
The LBG cB58 at z = 2.72 (Ellingson et al., 1996; Frayer et al., 1997; Seitz 
et al., 1998; Pettini et al., 2000), which is magnified strongly (by a factor of 
10-20) by a foreground cluster of galaxies at z = 0.37, and is at least ten 
times brighter than a typical LBG, was detected in the mm and submm by 
Baker et al. (2001) and van der Werf et al. (2002). However, after correcting 
for lensing, its 850-/im flux density is only about 0.1 mJy, below the level of 
confusion noise in SCUBA images, and similar to the flux density level of the 
statistical detection of high-redshift LBGs in the 850-/xm SCUBA image of 
the HDF-N (Peacock et al., 2000). In the field surrounding an overdensity 
of LBGs at z — 3.09, Chapman et al. (2001a) were successful in detecting 
bright submm emission that appears to be associated with diffuse sources of 
Lyman-a emission at the redshift of the overdensity, but were not included in 
the Lyman-break catalog. A key point to note is that the limits on submm- 
wave emission from LBGs are typically lower than expected if the relationship 
between UV spectral slope and far-IR luminosity observed for low-redshift 
low- luminosity starburst galaxies (Meurer et al., 1999) continues to high red- 
shifts. Goldader et al. (2002) indicate that the relationship does not appear 
to hold for the most luminous galaxies. 

The required sensitivity for successful submm observations of typical LBGs 
seems to be deeper than can be achieved using existing instruments. Obser- 
vations using future very sensitive, high-resolution interferometers certainly 
ALMA, and perhaps CARMA and SMA, will shed more light on the submm- 
LBG connection. 

The advent of the current generation of very sensitive X-ray observatories, 
Chandra and XMM-Newton, is generating a large sample of faint, hard X-ray 
sources, the luminosity of which is assumed to be dominated by high-redshift 
AGN (Fabian, 2000). Absorption and Compton scattering in large column den- 
sities of gas preferentially depletes soft X-rays, hardening the X-ray SEDs of 
gas-rich AGN. Such a population of hard, absorbed X-ray sources is required in 
order to account for the cosmic X-ray background radiation spectrum, which 
is harder than the typical SEDs of individual low-redshift AGN (Fabian and 
Barcons, 1992; Hasinger et al., 1996). Observations of the limited areas of 
the sky where both submm and X-ray data are available (Fabian et al., 1999; 
Hornschemeier et al., 2000; Mushotzky et al., 2001; Almaini et al., 2002) have 
tended to show little direct overlap between the X-ray and submm galaxies, 
although there are examples of X-ray-detected submm-wave galaxies (Bautz 
et al., 2000). The combined results of Bautz et al. and Fabian et al. reveal 
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that 2 out of 9 SCUBA galaxies are detected by Chandra. In the larger-area 
brighter 8-mJy survey, Almaini et al. (2001) identify only 1 out of 17 SCUBA 
galaxies using Chandra. Page et al. (2002) discuss further the submm prop- 
erties of X-ray sources. Perhaps of order 10% of known submm galaxies have 
faint hard X-ray counterparts that would be typical of dust-enshrouded AGN. 
There is also a statistical detection of excess submm-wave emission from the 
positions of faint high-redshift hard X-ray sources (Barger et al., 2001) and a 
positive submm-X-ray galaxy correlation function (Almaini et al., 2002). The 
lack of strong X-ray emission from a majority of submm galaxies lends cir- 
cumstantial support to the idea that much of their luminosity is derived from 
star formation and not from AGN accretion. However, some submm galaxies 
may have hydrogen column densities, and thus optical depths to Compton 
scattering, that are sufficiently great to obscure soft X-ray radiation entirely 
(> 10 24 cm~ 2 ). Even if they contained powerful AGN, these submm galaxies 
would be very faint in Chandra surveys, which reach detection limits of order 
10~ 17 ergcm~ 2 s^ 1 at soft 0.5-2keV wavelengths (Giacconi et al., 2002). They 
may be found in very deep observations using the the greater collecting area of 
XMM-Newton for hard X-ray photons. However, note that the 15-arcsec res- 
olution of XMM-Newton leads to confusion due to unresolved faint sources in 
the beam that is likely to impose a practical limit of order 10~ 15 erg cm -2 s _1 to 
the depth of a survey in the hard X-ray 2-8 keV band (Barcons et al., 2002). 
Deconvolution of the images from joint XMM-Newton/ Chandra deep fields, 
exploiting the sub-arcsec positional information from Chandra, will perhaps 
allow this limit to be exceeded. 

Finally, galaxies detected in far-IR surveys using ISO (for example Puget et 
al., 1999) out to redshifts z ~ 1 have been targeted for SCUBA submm obser- 
vations (Scott et al., 2000). The large arcmin-scale observing beam in 170-/xm 
ISO surveys makes identification of submm counterparts difficult, but progress 
has been made by combining sub-arcsec resolution deep radio images. The re- 
sults include some sources with apparently rather cool dust temperatures of 
order 30 K (Chapman et al., 2002d), and are generally consistent with red- 
shifts less than unity and dust temperatures of less than about 50 K for most 
of the galaxies. 



2.9 Alternative strategy for deep submm surveys 

An alternative strategy for finding submm galaxies has also been tried: tar- 
geted submm observations of faint radio-selected galaxies with accurate posi- 
tions and radio spectral index information detected at a 1.4-GHz flux density 
level brighter than 40-/iJy using the VLA (for example Richards, 2000). The 
radio source population at these faint flux density levels is expected to include 
mostly high-redshift star-forming galaxies, and only a minority of sources with 
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the more powerful synchrotron-emitting jets and lobes associated with parti- 
cles accelerated by AGN. By selecting those faint radio sources with faint 
K-band counterparts (Barger et al., 2000; Chapman et al., 2001b), it should 
be possible to sift out low-redshift galaxies from the target sample, and so 
generate a concentrated sample of luminous high-redshift galaxies to study in 
the submm. From Fig. 4, it is clear that a high-redshift faint 20-100-/zJy non- 
AGN 1.4-GHz radio source is likely to be a very luminous galaxy. Chapman 
et al. (2001b) claim that high-redshift submm galaxies can be detected using 
SCUBA at a rate of about one every hour using this method, which is signifi- 
cantly more rapid than the rate of about one every 10 h achieved in blank-field 
surveys. The efficiency of submm detection of optically-faint 1.4-GHz 20-/xJy 
radio galaxies at an 850-/im flux density greater than about 5 mJy is of order 
30-40%. Barger et al. (2000) detect 5/15 with / > 24 at 6mJy, and Chapman 
et al. (2001b) detect 20/47 with / > 25 at 4.5 mJy. 

The difficulty comes in the interpretation of the results. The additional selec- 
tion conditions of requiring first a radio detection, and then a faint optical 
counterpart will inevitably lead to the omission of a certain fraction of the 
submm galaxy population. Galaxies missing would include both the order 
of 15% of distant submm-selected galaxies that have relatively bright optical 
counterparts — the 'Class-2' SCUBA galaxies (Ivison et al., 2000a; Smail et al., 
2002) — and any very high-redshift submm galaxies that cannot be detected 
at the VLA, despite lying on the far-IR-radio correlation (Condon, 1992; see 
Fig. 8). There is also likely to be a bias (at all redshifts) towards both detecting 
AGNs, in which radio flux densities are boosted above the level expected from 
the standard radio-far-IR correlation, and the exclusion of a small population 
of low-redshift submm galaxies, which would be too bright in the radio and 
-fTband to be included in the survey. The rare, and perhaps especially inter- 
esting, submm galaxies at the lowest and highest redshifts are thus likely to 
be missing from radio pre-selected surveys, as are the distant submm galaxies 
with brighter optical magnitudes that are likely to be easiest to follow-up and 
investigate. 

The size of these selection effects is difficult to quantify at present. However, 
their existence can be inferred from the diverse multiwaveband properties dis- 
played by submm-selected galaxies from imaging surveys, a significant fraction 
of which are known not to be detected in deep radio observations (Downes et 
al., 1999; Smail et al., 2000). In some cases, the depth of the radio images used 
to make the comparison could be improved by a factor of several; it is possi- 
ble, but we suspect unlikely, that all submm sources have radio counterparts 
lurking just below existing detection thresholds. 

It is instructive to compare the results of the blank-field and faint radio pre- 
selected submm surveys. In a true 260-arcmin 2 blank-field survey, Scott et al. 
(2002) found surface densities of 550^0 an d 180 ± 60 deg -2 submm galaxies 
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brighter than 6 and lOmJy respectively, while Borys et al. (2002) estimate 
164 ± 28 cleg -2 brighter than 12mJy in a 125-arcmin 2 survey, with a rather 
conservative error estimate from 12 detections. The surface densities resulting 
from the faint radio-selected investigations of Barger et al. (2000) and Chap- 
man et al. (2001b) are 430 deg -2 brighter than 4mJy and 135 deg -2 brighter 
than 10 mJy.[^] This indicates incompleteness in the faint radio-selected counts 
by at least about 25%. This is perfectly acceptable for gathering a list of high- 
redshift galaxies for further study; however, for statistical analysis of a large 
sample of many tens of submm galaxies, the uncertainty introduced in the de- 
rived properties of the population due to incompleteness is expected to domi- 
nate the Poisson uncertainty, and so limit the accuracy of the inferences that 
can be made. Once the bright counts of distant 850-/iin galaxies are known 
from direct, hopefully unbiased, large-area SCUBA imaging surveys (Fox et 
al., 2002; Scott et al., 2002) and matched with very deep radio images (Ivison 
et al., 2002), then the value of this shortcut for compiling a large sample of 
submm galaxies can be assessed. The same strategy would be very valuable for 
wide-field MAMBO 1.2-mm surveys (Dannerbauer et al., 2002). At present, 
some care needs to be taken in the interpretation of faint radio pre-selected 
surveys. 



2.10 Determining redshifts of submm galaxies 

By virtue of their uniform selection function with redshift, it is impossible to 
even indicate the redshift of submm galaxies from single- wavelength submm 
flux densities alone. The relatively poor positional accuracy of submm images 
also makes it very difficult to identify an unambiguous optical counterpart 
for spectroscopic follow-up. This will be discussed in more detail when the 
properties of individual galaxies are described in Section 3. Here we discuss 
some general features of the prospects for determining their redshifts from 
submm, far-IR and radio observations. 

2.10.1 Photometric redshifts from far-IR SEDs 

The submm-far-IR SED of dusty galaxies is thermal, and so redshifting a 
fixed template SED affects the observed colors in exactly the same way as 
changing the dust temperature. This means that even when multi-frequency 
far-IR data is available, the redshift will be uncertain, unless there is informa- 
tion about the intrinsic dust temperature of the source. This effect is related 



11 Barger et al., (2000) also impose a limit of less than 840 deg~ 2 brighter than 
6mJy at 850-//m, based on serendipitous detections of sources in their 2.5-arcmin- 
wide SCUBA images around the AGN radio galaxies. 
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to the selection effect in favor of cold galaxies illustrated in Fig. 5. The ex- 
pected colors of dusty galaxies, as a function of (1 + z)/T d , the parameter 
that can be constrained in light of this degeneracy, is shown in Fig. 6. Colors 
in a variety of observing bands in both submm atmospheric windows and for 
the observing bands of the MIPS instrument on the Space InfraRed Telescope 
Facility (SIRTFj^j satellite are included. It is vital to stress that without 
knowledge of the dust temperature, it is impossible to determine a redshift 
from any combination of multicolor broadband far-IR/submm data. The de- 
generacy between T d and z is lifted partially by combining information derived 
from radio observations, if T d > 60 K. It may also become clear that there is 
a Universal temperature-luminosity relation that extends to high redshifts, 
and can be exploited to determine redshift information using a far-IR-submm 
color-magnitude relation. In the absence of such a Universal relation, far-IR 
and submm colors can only be used to fix the parameter (1 + z)/T d reliably. 

2.10.2 Radio- submm photometric redshifts 

In Fig. 4, an estimate of the flux density of radio emission of a template dusty 
galaxy was shown as a function of redshift. We now investigate the radio 
properties of submm galaxies, which provide useful information about their 
SED and redshift. There is an excellent observed correlation between the radio 
and far-IR (60- and 100-/im) flux densities of low-redshift galaxies over 4 
orders of magnitude in luminosity, reviewed by Condon (1992) and recently 
investigated out to z ~ 0.3 by Yun et al. (2001). If this correlation is assumed 
to hold to high redshifts, then submm-selected galaxies should be detectable in 
the deepest 10-/zJy-RMS 1.4-GHz VLA radio images out to redshifts of order 
3. Note that as a result of this correlation, optical spectroscopy of faint non- 
AGN radio galaxies alone can be used to trace the evolution of star-formation 
activity to redshifts z ~ 1.2, beyond which spectroscopic redshifts are hard to 
determine (Haarsma et al., 2000). 

The far-IR-radio correlation is thought to be due to a match between the 
rate at which the optical/UV radiation from young stars is absorbed by dust 
on a local scale in star-forming regions of galaxies, and re-emitted as ther- 
mal far-IR radiation, and the radio luminosity from the same regions (Harwit 
and Pacini, 1975). The radio luminosity is due to both free-free emission in 
Hn regions, and more importantly at frequencies less than about 10 GHz, to 
the synchrotron emission from relativistic electrons accelerated in supernova 
shocks. If an AGN is present, then it is likely that its accretion disk will 
provide an additional source of UV photons to heat dust, and both the disk 
and outflows will generate shocks to accelerate relativistic electrons. There is 
little reason to expect these effects to be proportional, unlike UV heating by 



sirtf.caltech.edu 



30 



o 

• H 
-+-> 

• i — i 

m 

G 

0) 
T3 

X 



Black lines: 1.4- GHz 
at 10, 20, 40, 80, 160K 




100 



10 



0.1 



0.01 



10" 



10 



-4 



0.5 1 2 



10 



(l + z)(7 1 /40K) 



-l 



Fig. 6. The ratio of flux densities expected in different observing bands as a function 
of the degenerate redshift /dust temperature parameter, compared with the flux 
density expected in the 70-/mi band at which the SIRTF satellite will be very 
sensitive. Where the lines have steep gradients, measured colors from multi-band 
data locate the peak of the dust SED accurately in the observer's frame, providing 
a measurement of temperature-redshift. The degeneracy between and z can be 
lifted slightly by including radio data (Blain, 1999a; Yun and Carilli, 2002), if the 
dust temperature is greater than about 60 K (see also Fig. 7). If deep near-IR and 
optical images can be included, and the optical counterpart to the galaxy can be 
readily identified, then conventional photometric redshifts can be determined from 
stellar synthesis models. However, care must be taken as it is unclear whether the 
SEDs of very dusty galaxies have familiar restframe-optical spectral breaks. 

massive stars and particle acceleration by supernova shocks in star-forming re- 
gions. Radio-quiet QSOs tend to lie on the radio-loud side of the far-IR-radio 
correlation, while radio-loud AGN can lie up to three orders of magnitude 
away. 



Carilli and Yun (1999, 2000) demonstrated that the radio-submm color is a 
useful redshift indicator, assuming that dusty galaxies have simple synchrotron 
SEDs in the radio waveband and thermal dust spectra in the submm and far-IR 
wavebands. The radio-submm color was also considered in the interpretation 
of the redshifts of galaxies detected in the submm surveys by Hughes et al. 
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Fig. 7. The behavior of the Carilli-Yun 1.4-GHz to 850- /^m radio-submm redshift 
indicator. The left panel shows the ratios of 1.4-GHz:850-^m flux density predicted 
from empirical SEDs by Carilli and Yun (1999; dotted lines) and Dunne et al. (2000; 
dashed line; see also Carilli and Yun, 2000). Predictions for the ratio based on the 
results of Blain (1999a) are also shown assuming radio-far-IR SEDs with various 
dust temperatures, but which all lie on the far-IR-radio correlation (Yun et al., 2001; 
solid lines). The flux ratio is a good indicator of redshift, clearly separating high- and 
low-redshift galaxies. Both synchrotron and free-free radio emission are included, 
and the dust temperature and radio properties evolve with redshift self-consistently, 
as modified by the CMB. In the right panel, the solid curves are replotted as a 
combined function of temperature and redshift, emphasizing that for Td < 60 K, the 
inferred temperature and redshift are degenerate, just as for a thermal spectrum 
(Fig. 6). For Td > 60 K the flux ratio becomes a non-degenerate redshift indicator. 

(1998), Lilly et al. (1999) and Eales et al. (2000). The Carilli-Yun redshift 
indicator is subject to a degeneracy between dust temperature and redshift, 
at dust temperatures less than about 60 K: see Fig. 7 in which its form is shown 
for a variety of SEDs. Hot, distant galaxies are difficult to distinguish from 
cool, low-redshift ones (Blain, 1999a, b). Despite this degeneracy, the Carilli- 
Yun redshift indicator is very useful, especially for investigating optically faint 
submm galaxies for which almost no other information is available (Smail 
et al., 2000). AGN, which are expected to be radio- loud as compared with 
the standard far-IR-radio correlation, lead to conservatively low Carilli-Yun 
estimated redshifts. 

At the highest redshifts, two additional factors arise to modify the relation 
expected. First, synchrotron emission from relativistic electrons is likely to be 
suppressed. The intensity of synchrotron emission depends on only the en- 
ergy density in the interstellar magnetic field, while the total cooling rate of 
the electrons depends on the sum of the energy densities in the interstellar 
magnetic field, for synchrotron losses, and in the ISRF, for inverse Compton 
scattering losses. The tight low-redshift far-IR-radio correlation implies that 
these energy densities are proportional over a very wide range of galaxy prop- 
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erties (Volk, 1988). However, above some critical redshift, the energy density 
in the CMB will always rise to dominate the ISRF, upsetting this balance. 
Inverse- Compton electron cooling will then dominate and the amount of syn- 
chrotron emission will decline. Note, however, that free-free emission does not 
suffer this suppression at high redshifts, and so should remain detectable out 
to any practical redshift. The almost flat SED of free-free emission ensures 
a more favorable K correction for high-redshift galaxies than expected for a 
pure synchrotron emission spectrum. The free-free emission spectrum cuts off 
only at photon energies greater than the thermal energy of emitting electrons 
in Hn regions. In the absence of free-free optical depth effects, for 10 4 " 5 K this 
corresponds to optical frequencies (Yun and Carilli, 2002). In Fig. 8 we show 
the effects of CMB suppression of synchrotron emission, for a ratio of energy 
densities in the magnetic field and the ISRF of 0.33, which is reasonable for 
M82 and the Milky Way (Hummell, 1986). We assume a galaxy SED template 
that lies on the standard far-IR-radio correlation at z = (Condon, 1992). 

Secondly, again most significant at high redshifts, a minimum dust tempera- 
ture is imposed by the rising CMB temperature: dust must be hotter than the 
CMB. Given that observed dust temperatures in ultraluminous galaxies seem 
to lie in the range 40-100 K, then this may become an important factor at 
redshifts 10 < z < 30, if an early generation of stars generates the heavy ele- 
ments required to form dust prior to these redshifts. For cooler Milky- Way- like 
SEDs, this effect would be important at z ~ 5. An increase in dust temper- 
ature due to CMB heating at high redshifts shifts the peak of the dust SED 
to higher frequencies, counteracting the beneficial K correction illustrated in 
Fig. 4. As a result, there is a firm upper limit to the redshift at which submm 
continuum radiation can be exploited to image the most distant galaxies effi- 
ciently, even if these galaxies do contain dust. This is illustrated by the flux 
density-redshift relations at 230, 90 and 30 GHz shown in Fig. 8: at redshifts 
greater than about 15 all three curves have the same form. 

Both of these factors were included in the derivation of the curves in Fig. 7, 
which should thus provide an accurate guide to the usefulness of the Carilli- 
Yun redshift indicator out to the highest redshifts. Note that the indicator 
becomes of little use for the most distant galaxies, for which an almost constant 
radio:submm flux ratio of about 10~ 3 is expected. In addition, the hottest 
dusty galaxies may be more likely to contain AGN, and thus to lie on the 
radio-loud side of the far-IR-radio correlation. This effect could make the 
interpretation of Fig. 7 for determining redshifts ambiguous, even if > 60 K. 
However, it should always guarantee a conservative estimate of the redshift 
for any observed galaxy. 

Foreground absorption is not a problem for very high-redshift sources of any 
radio, submm or far-IR radiation. It is thus likely that the most sensitive 
future instruments at submm and radio wavelengths, ALMA and the SKA 
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Fig. 8. Some key features of flux density-redshift relations expected at a range 
of wavelengths, extending to very high redshifts. CMB heating of dust at z > 10 
prevents the mm- wave K correction from assisting the detection of very high redshift 
galaxies: the flux density-redshift relation has the same redshift dependence beyond 
z ~ 15 at 230, 90 and 30 GHz. CMB cooling of relativistic electrons suppresses 
synchrotron radio emission beyond z ~ 5, as shown by the thin dashed line. Realistic 
free-free emission is also included in the model represented by the thick dashed line 
(Condon, 1992; Yun et al., 2001), significantly increasing the radio emission expected 
from very high redshifts. An estimate of the flux density from a 3 x 10 4 -K stellar 
photosphere at 5/^m is also shown, cutoff at the redshift beyond which Lyman-a 
absorption is redshifted through the band. Note that there is probably a maximum 
redshift above which dust does not exist, and so beyond which thermal emission 
from dust can never be detected; this effect is not included here. 

will both be able to detect 'first light' galaxies. Note, however, that a mask of 
foreground structure may be significant for radio observations (Waxman and 
Loeb, 2000). A practical limit to the capability of determining the history of 
very early star formation from an 'SKA deep field' could also be set by the low, 
cosmologically-dimmed surface brightness of galaxies at the highest redshifts, 
and their potentially overlapping emission regions. The importance of both of 
these factors is expected to depend critically on the unknown physical sizes 
of the first galaxies. ALMA will probably be limited to observe galaxies out 
to a maximum redshift set by the requirement that sufficient metals have 
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been generated to form obscuring dust at that early epoch. To compare what 
might be possible in the near-IR waveband, the emission from a 3 x 10 4 -K 
blackbody stellar photosphere at 5 //m is shown in Fig. 8, cutoff at the redshift 
beyond which absorption by redshifted Lyman-a becomes important in the 
band. This indicates the potential for probing the earliest galaxies using a 
near-IR camera on the Next Generation Space Telescope (NGST). Free-free 
emission and redshifted near-IR stellar emission may thus be the best routes 
to the detection of the first galaxies at redshifts greater than 10, if they exist. 



3 The observed properties of submm-selected galaxies 

Well over 100 submm-selected galaxies are now known (see Fig. 9), although 
their redshifts and detailed astrophysical properties are very largely uncertain. 
The key information available about their properties comes from observations 
of discrete galaxies made using the SCUBA and MAMBO bolometer array 
cameras at wavelengths of 450, 850 and 1200 /im. Counts of distant galaxies 
at far-IR wavelengths of 95 and 175 fim have also been measured using the 
PHOT instrument aboard ISO. Limits to the counts at 2.8 mm have been 
obtained using the Berkeley-Illinois-Maryland Association (BIMA) mm-wave 
interferometer. The results of all the relevant observations are compiled in 
Figs. 9 and 10. Information is also available about the population of mid-IR 
15-/im sources using the CAM instrument aboard ISO (Altieri et al., 1999; 
Elbaz et al., 1999): see Fig. 22. 

Measurements of the observed intensity of background radiation from the radio 
to far-UV wavebands are shown in Fig. 11. The background measurements 
include lower limits obtained by summing over the observed counts plotted 
in Figs. 9 and 10. the deepest counts at 850 and 15 fim come from surveys 
made in fields magnified by gravitational lensing clusters. Surface brightness 
conservation in the lensing process ensures that the mean background intensity 
in the direction of a lens should be the same as that in a blank field. 

From the properties of the counts and backgrounds alone, without any details 
of the individual galaxies involved, it is possible to infer important details 
about the population of distant dust-enshrouded galaxies. 

The significant surface density of the faint SCUBA and MAMBO galaxies, 
when coupled to plausible SEDs (Blain et al., 1999b; Trentham et al., 1999; 
Dunne et al., 2000), clearly indicates that the luminosity function of dis- 
tant dusty submm galaxies is much greater than that of low-redshift IRAS 
galaxies (Saunders et al., 1990; Soifer and Neugebauer, 1991), and undergoes 
very strong evolution. An extrapolation of the low-redshift luminosity function 
without evolution predicts a surface density of galaxies brighter than 5mJy 
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Fig. 9. A summary of count data from several mm, submm and far-IR surveys. The 
overplotted curves are derived in models that provide good fits to the compilation 
of data, and are updated from the results in the listed MNRAS papers (Blain et al., 
1999b, c). Identical symbols represent post-1999 data from the same source. The 
errors are shown as la values unless stated. The 2.8 mm data (square) is from Wilner 
and Wright (1997). At 850 /im, in order of increasing flux (less than 15 mJy), data is 
from Blain et al. (1999a), Cowie et al. (2002) with 90% confidence limits; Hughes et 
al. (1998), Chapman et al. (2001b), Barger et al. (1999a, 1998), Smail et al. (1997), 
Eales et al. (1999, 2000) consistent with the increased area reported by Webb et 
al. (2002a), Borys et al. (2001), Barger et al. (2000) and Scott et al. (2002). The 
data points between about 2 and 10 mJy are consistent with a steep integral source 
count iV(> S) oc S a , with a power-law index a ~ —1.6. The counts at brighter flux 
densities are likely to steepen considerably; note that the counts must turn over at 
fainter flux densities to have a < — 1 to avoid the background radiation intensity 
diverging. At 95 /im, the data is from Juvela et al. (2001), Kawara et al. (1998), 
Matsuhara et al. (2000), Serjeant et al. (2001) and Linden- Vornle et al. (2000). 

at 850 /iin of only about 0.25 deg -2 , as compared with the observed density 
of several 100 deg -2 (Fig. 9). Because of the flat flux density-redshift relation 
in the submm shown in Fig. 4, a 5-mJy SCUBA galaxy at any moderate or 
high redshift (z > 0.5) has a luminosity greater than about 8 x 10 12 L . Imme- 
diately, this tells us that the comoving density of high-redshift galaxies with 
luminosities in excess of about 10 13 L is 400 times greater than at z — 0. 
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Fig. 10. Counterpart to Fig. 9 for three other observing bands. The data at 1.2 mm 
(circles at flux densities less than lOmJy) are from Bertoldi et al. (2001), Carilli et 
al. (2000) and Carilli (2001). The data at 450 /im (circles at 10-50 mJy) are from 
Smail et al. (2002), with limits from Smail et al. (1997) and Barger et al. (1998). 
The data at 175 p ( 100 mJy) are from Kawara et al. (1998), Puget et al. (1999), 
Matsuhara et al. (2001), Juvela et al. (2001), Dole et al. (2001) and Stickel et al. 
(1998). 

We stress that the submm K correction ensures that the redshift has little 
effect on the results: the count would be approximately the same whether the 
population is concentrated at z ~ 1 or extends from z ~ 2 to 10. 

This estimate is subject only to an uncertainty in the dust temperature, which 
is assumed to be about 40 K. Even if the dust temperature of some of the 
galaxies is as low as the 20 K found for low-redshift spiral galaxies, then their 
luminosity is still about 8 x 1O 11 L , considerably greater than the several 
10 10 L expected for typical spiral galaxies. This issue can be addressed by 
taking into account both the observed background spectrum and the counts 
at different wavelengths. 

The submm-wave background radiation spectrum can also be exploited to 
provide information about the form of evolution of the luminosity function. 
The submm-wave background, measured directly using COBE-FIRAS (Puget 
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Fig. 11. The observed intensity of cosmic background radiation between the radio 
and far-UV wavebands. The great majority of the background energy density in 
the Universe derived from sources other than the CMB is represented in this figure. 
Almost all of the rest appears in the X-ray waveband. Some significant uncertainty 
remains, but the combination of measurements and limits indicates that a compa- 
rable amount of energy is incorporated in the far-IR background, which peaks at a 
wavelength of about 200 fxm, and in the near-IR/optical background, which peaks 
at a wavelength between 1 and 2 fim. The data originates from a wide range of 
sources: 1. Fixsen et al. (1998); 2. Puget et al. (1996); 3. Blain et al. (1999a); 4. 
Schlegel et al. (1998); 5. Hauser et al. (1998); 6. Lagache et al. (2000a) see also 
Kiss et al. (2001); 7. Puget et al. (1999); 8. Kawara et al. (1998); 9. Finkbeiner et 
al. (2001); 10. Stanev and Franceschini (1998); 11. Altieri et al. (1999); 12. Dwek 
and Arendt (1998); 13. Wright and Johnson (2002); 14. Pozzetti et al. (1998); 15. 
Bernstein (1999) and Bernstein et al. (2002); 16. Toller et al. (1987); 17. Armand et 
al. (1994); 18. Lampton et al. (1990); and 19. Murthy et al. (1999). For a detailed 
review of cosmic IR backgrounds see Hauser and Dwek (2001). Note that Lagache et 
al. (2000a,b) claim that the Finkbeiner et al. points (9) could be affected by diffuse 
zodiacal emission. Where multiple results are available in the literature the most 
sensitive result is quoted. 

et al., 1996; Hauser et al., 1999; Schlegel et al., 1999), reasonably exceeds the 
sum of the measured flux densities of discrete galaxies detected in SCUBA 
surveys (Smail et al., 1997, 2002; Blain et al., 1999a). However, the submm 
background makes up only a small fraction of the total energy density in 
the far-IR background, which peaks at a wavelength of about 200 /im and is 
generated by galaxies at redshift z ~ 1. The relatively flat source SEDs and 
the rate of change of the cosmic volume element at this redshift conspire to 
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generate most of the background light, just as in the radio, X-ray, optical 
and near-IR wavebands. The mm and submm background radiation is unique 
in originating at a higher redshift. Very little of the background is expected 
to be generated at z < 1, and so it is an important signature of high-redshift 
galaxy formation. Despite representing only a small fraction of the total energy 
density in the cosmic background radiation, the mm-wave background is one 
of the cleanest measures of activity in the distant Universe. 

There are significant consequences for the evolution of galaxies at high red- 
shifts due to the observed smooth power-law form of the background spec- 
trum, I v oc z/ 264 , for v < 500 GHz (Fixsen et al., 1998), which originates at 
moderate to high redshifts, on account of the submm-wave K correction. The 
shape of the background radiation spectrum at frequencies greater than about 
100 GHz can be approximated quite accurately by associating an evolving co- 
moving volume emissivity (pl in units of Wm -3 ) with an SED that peaks at 
a single frequency u , so that e u oc p^{z)5{u — z/ ) (Blain and Longair, 1993b), 
and then integrating over cosmic volume over a fixed angle on the sky. If the 
SED, via u Q , is assumed not to evolve strongly with redshift — there is no clear 
evidence that it does — then in order to reproduce the observed slope of the 
mm/submm background spectrum, Pl(z) oc (l + z)~~ 1A is required for z ^> 1, 
and so the comoving luminosity density of dust-enshrouded galaxies must de- 
cline at large redshifts. If it did not decline, then the background spectrum 
measured by COBE would be too flat, with too much energy appearing at long 
wavelengths. This argument has been made using Monte-Carlo simulations of 
Pl(z) by Gispert et al. (2000). A similar set of simulations have been carried 
out by Eales et al. (2000), taking into account the observed background ra- 
diation spectrum, counts and inferred redshift distribution of submm-selected 
galaxies. 

An approximately equal fraction of the cosmic background radiation energy 
density emerges in the near-IR/optical and far-IR wavebands (Fig. 11). Be- 
cause dusty galaxies do not dominate the total volume emissivity at low red- 
shifts (Sanders, 1999; Yun et al., 2001), then the volume emissivity of dusty 
galaxies must increase by a factor of at least 10, matching the significant 
evolution of the population of galaxies observed in the optical waveband at 
z < 1 (Lilly et al., 1996), to avoid the intensity of the far-IR background 
radiation being significantly less than observed. Only a very small fraction of 
the total far-IR luminosity from all low-redshift galaxies comes from galax- 
ies more luminous than 1O 12 L , yet as discussed above in the context of the 
submm-wave counts, these luminous galaxies are much more numerous at high 
redshifts, by a factor of several hundred. These twin constraints demand that 
the form of evolution of the luminosity function of dusty galaxies cannot be 
pure density evolution, a simple increase in the comoving space density of all 
far-IR-luminous galaxies. If the counts were to be reproduced correctly in such 
a model, then the associated background radiation spectrum would be much 
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greater than observed. A form of evolution similar to pure luminosity evolu- 
tion, in which the comoving space density of galaxies remains constant, but 
the value of L*, the luminosity that corresponds to the knee in the luminosity 
function, increases — in this case by a factor of order 20 — is consistent with 
both the submm-wave counts and background intensity. 

By a more rigorous process, taking into account all available information, 
including the need to normalize the results to the observed low-redshift pop- 
ulation of dust-enshrouded galaxies from the IRAS luminosity function and 
the populations of galaxies observed by ISO at z ~ 1, the evolution of the 
luminosity density pl can be constrained. The results have been discussed by 
Blain et al. (1999b, c), as updated in Smail et al. (2002), and by Eales et al. 
(2000). They are discussed further in Section 5 below. 



3. 1 Confusion 

Source confusion, the contribution to noise in an image due to the super- 
imposed signals from faint unresolved sources clustering on the scale of the 
observing beam (Condon, 1974; Scheuer, 1974), is a significant problem for 
observations in the submm waveband (Blain et al., 1998; Eales et al., 2000; 
Hogg, 2001). This is due to the relatively coarse (~ 10 arcsec) spatial resolution 
currently available. In fact, a significant fraction of the noise in the deepest 
850-/im SCUBA image of the HDF-N (Hughes et al., 1998) can be attributed 
to confusion (Peacock et al., 2000). At present, the practical confusion limit 
for galaxy detection in SCUBA observations at the atmospherically favored 
850-/im wavelength is about 2 mJy This limit makes it difficult to determine 
accurate sub-arcsec positions for the centroids of the submm emission from 
faint SCUBA-selected galaxies, rendering follow-up observations more chal- 
lenging. Unfortunately, experience has shown that many known high-redshift 
galaxies, especially optically-selected LBGs, are typically fainter than the con- 
fusion limit, and so are difficult to study using SCUBA. 

The variety of count data for dusty galaxies shown in Figs. 9 and 10 can be 
used to estimate the effect of source confusion in observations made at a wide 
range of frequencies and angular scales. 

The distribution of flux density values from pixel to pixel in an image due 
to confusion noise depends on the underlying counts of detected galaxies. 
Confusion noise is always an important factor when the surface density of 
sources exceeds about 0.03 beam -1 (Condon, 1974). The results of a confusion 
simulation for the 14-arcsec SCUBA 850-/im beam are shown in Fig. 12: see 
also Hogg (2001), Eales et al. (2000) and Scott et al. (2002). Note that Eales et 
al. assume a very steep count and obtain larger values of confusion noise than 
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Fig. 12. Histograms showing the simulated effects of confusion noise in deep SCUBA 
integrations at 850 jum. Left: the expected distribution of pixel flux densities when 
the telescope samples the sky in a standard (—0.5,1,-0.5) chopping scheme, with no 
additional noise terms present. The flux distribution is non-Gaussian, with enhanced 
high- and low-flux tails as compared with the overplotted Gaussian, which has 
the width predicted by simple calculations (Fig. 13). Right: the same confusion 
noise distribution is shown in the right-hand panel, but convolved with Gaussian 
instrument and sky noise with an RMS value of 1.7 mJy, which is typical of the 
noise level in the SCUBA Lens Survey (Smail et al., 2002). At this noise level, the 
additional effect of confusion noise is small. 

those shown here. Observations made with finer beams at the same frequency 
suffer reduced confusion noise, while for those made in coarser beams the 
effects are more severe: compare the results for the much larger 5-arcmin beam 
in the three highest frequency submm-wave observing bands of the planned 
Planck Surveyor space mission all-sky survey shown in Blain (2001a). 

The simulated confusion noise distribution is non-Gaussian (Fig. 12), but can 
be represented quite accurately by a log-normal distribution, leading to many 
more high-flux-density peaks in an image than expected assuming a Gaussian 
distribution of the same width. The width of the central peak of the distri- 
bution in flux density is approximately the same as the flux density at which 
the count of sources exceeds 1 beam" 1 . This provides a useful indication of the 
angular scales and frequencies for which confusion noise is likely to be signifi- 
cant, and of the limit imposed to the effective depth of surveys by confusion 
for specific instruments: see Fig. 13. 

3.1.1 Confusion and follow-up observations of submm galaxies 

The real problem of confusion for identifying and conducting multiwaveband 
studies of submm-selected galaxies is illustrated by the results of the first gen- 
eration of surveys. The very deepest optical image that matches a submm-wave 
survey is the HDF-N, in which there are several tens of faint optical galaxies 
(at R > 26) that could be the counterpart to each SCUBA detection (Hughes 
et al., 1998; Downes et al., 1999a). It is thus impossible to be certain that 
a correct identification has been made from the submm detection image and 
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Fig. 13. An approximate measure of the 1-a confusion noise expected as a function 
of both observing frequency and angular scale from the mm to mid-IR waveband 
(updated from Blain et al., 1998). The contributions from extragalactic and Galac- 
tic sources are shown in the left and right panels respectively. Radio-loud AGN 
may make a significant contribution to the top left of the jagged solid line (Tof- 
folatti et al., 1998). A Galactic cirrus surface brightness of Bq = lMJysr^ 1 at 
100 fim is assumed. The ISM confusion noise is expected to scale as Bq 5 (Helou 
and Beichman, 1990; Kiss et al., 2001). The bands and beamsizes of existing and 
future experiments (see Tables 1 and 2) are shown by: circles — Planck Surveyor; 
squares — BOOMERANG; empty stars — the SuZIE mm-wave Sunyaev-Zeldovich 
instrument; triangles — BOLOCAM, as fitted to CSO (upper 3 points) and the 
50-m Large Millimeter Telescope (LMT; lower 3 points); filled stars — SCUBA 
(and SCUBA-II); diamonds — Herschel; asterixes — Stratospheric Observatory for 
Infrared Astronomy (SOFIA); crosses — SIRTF. The resolution limits of the inter- 
ferometric experiments ALMA and SPECS lie far below the bottom of the panels. 
The confusion performance of the 2.5-m aperture BLAST balloon-borne instru- 
ment is similar to that of SOFIA. Confusion from extragalactic sources is expected 
to dominate over that from the Milky Way ISM for almost all of these instruments. 



optical data alone: compare the identifications in Smail et al. (1998a, 2002). In 
some cases, EROs and faint non-AGN radio galaxies (Smail et al., 1999, 2000; 
Gear et al., 2000; Lutz et al., 2001) can be associated with submm galaxies, 
especially after higher-resolution mm-wave interferometry observations have 
provided more accurate astrometry for the submm detection (Downes et al., 
1999; Frayer et al., 2000; Gear et al, 2000; Lutz et al, 2001), to reduce the 
effects of submm confusion, with the investment of significant amounts of ob- 
serving time. The surface density of both EROs and faint non-AGN radio 
galaxies is less than that of the faintest optical galaxies, and so the probabil- 
ity of a chance coincidence between one and a submm galaxy is reduced. A 
very red color and detectable radio emission from a high-redshift galaxy are 
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both likely to indicate significant star- formation/ AGN activity and/or dust 
extinction, making such galaxies better candidate counterparts even in the 
presence of confusion-induced positional uncertainties. 



3.2 Multi-waveband follow-up studies 

A great deal of telescope time has been spent so far to detect and study 
submm-selected galaxies in other wavebands. In many cases, rich archival data 
predated the submm observations: most notably in HDF-N (Hughes et al., 
1998). Considerable data was also available in the fields of rich clusters (Smail 
et al, 1997, 1998a, 2002; Cowie et al, 2002), in the region of the Eales et al. 
(1999, 2000) surveys in Canada-France Redshift Survey (CFRS) fields, which 
include the Groth Strip, and in the deep Hawaii survey fields (Barger et al., 
1999a). The results of follow-up deep optical and near-IR imaging (Frayer et 
al., 2000) and spectroscopy (Barger et al., 1999b), mm- wave continuum imag- 
ing (Downes et al., 1999a; Bertoldi et al., (2000); Frayer et al. 2000; Gear et 
al. 2000; Lutz et al., 2001; Dannerbauer et al., 2002) and molecular line spec- 
troscopy (Frayer et al., 1998, 1999; Kneib et al., 2002) have been published, 
and many additional studies are under way. The time spent following up the 
SCUBA Lens Survey (Smail et al., 2002) exceeds by almost an order of mag- 
nitude the time required to make the submm discovery observations (Smail 
et al., 1997). The difficulty of the task is highlighted by the identification of 
plausible counterparts to these galaxies being only about 60% complete over 
4y later at the start of 2002. The follow-up results from a well-studied subset 
of galaxies in the SCUBA Lens Survey are shown in Figs. 14-20, in order of de- 
creasing 850-/im flux density. These are chosen neither to be a representative 
sample of submm galaxies, nor to be a sufficiently large sample for statistical 
studies, but rather to present a flavor of the range of galaxies that can be de- 
tected in submm-wave surveys for which high-quality multi-waveband data is 
available. The galaxies that are presented are chosen to have good positional 
information, and redshifts where possible. Other detections for which excellent 
mult i- wavelength follow-up data are available include the brightest source in 
the HDF-N SCUBA image (Hughes et al., 1998; Downes et al., 1999a), an 
ERO detected in the CUDSS survey by Eales et al. (1999) (Gear et al., 2000), 
a z = 2.8 QSO in a cluster field (Kraiberg-Knudsen et al., 2001), and the 
substantially overlapping catalogs of galaxies detected by Cowie et al. (2002) 
in deeper images of 3 of the 7 clusters in the Smail et al. lens survey. 

3. 2. 1 Optical /near-IR 

The properties of submm galaxies in the optical waveband, corresponding to 
the rest-frame UV waveband, appear to be very diverse (Ivison et al., 2000a). 
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This may be due in part to their expected broad redshift distribution. How- 
ever, given that two submm galaxies at the reasonably high redshifts z = 2.5 
and 2.8 are known to be readily detectable at B ~ 23 (Ivison et al., 1998, 2000, 
2001) — before correcting each for the magnitude (factor of about 2.5) of ampli- 
fication due to the foreground cluster lens — while most others are very much 
fainter (Smail et al., 2002; Dannerbauer et al., 2002), it is likely that much 
of the spread in their observed properties is intrinsic. As most counterparts 
are extremely faint, confirmation of their nature requires a large, completely 
identified sample of submm galaxies with known redshifts, which is likely to be 
some time away. It is possible that optically faint submm galaxies have similar 
properties, an issue that can be addressed when deep near-IR observations are 
available. 

Ivison et al. (2000a) and Smail et al. (2002) proposed a 3-tier classification sys- 
tem to stress the varied nature of submm galaxies (three being an eminently 
sensible number of classes for 15 galaxies!). Class-0 galaxies are extremely 
faint in both the observed optical and near-IR wavebands. Class-1 galaxies 
are EROs, very faint in the optical but detectable in the near-IR, while Class- 
2 galaxies are relatively bright in both bands. It is unclear how closely this 
scheme reflects the underlying astrophysics of the submm galaxies; however, 
the classification separates the optically bright galaxies (Class 2s), for which 
the acquisition of optical redshifts and confirming CO redshifts are likely to 
be practical, and the fainter galaxies, for which this will be a great challenge 
(Class 0s). A similar approach for MAMBO sources has been discussed by 
Dannerbauer et al. (2002). Note that submm galaxies could change classifica- 
tion by having different redshifts, despite identical intrinsic SEDs. 



3.2.2 Ultradeep radio images 

The surface density of the faintest radio sources that can be detected using the 
VLA (Richards, 2000) is significantly less than that of optical galaxies, and 
so an incorrect radio counterpart to a submm-selected galaxy is relatively un- 
likely to be assigned by chance. If 1.4-GHz VLA images are available at a flux 
limit approximately 1000 times deeper than 850-/xm images, then the radio 
counterparts to non-AGN submm galaxies should be detectable to any red- 
shift (see Fig. 7). Deep radio follow-up observations of submm-selected galaxies 
yielding cross identifications have been discussed by Smail et al. (2000), and 
further information about very faint radio sources in the field of the UK 8-mJy 
SCUBA survey (Fox et al., 2002; Scott et al., 2002) should soon be available in 
Ivison et al. (2002). Despite an extremely deep radio image (Richards, 2000), 
the brightest submm galaxy detected in HDF-N (Hughes et al., 1998) does not 
have a radio detection, probably indicating a very high redshift. The survey 
results reported by Eales et al. (2000) and Webb et al. (2002a) were discussed 
in the context of radio data covering the same fields; however, the radio survey 
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is not deep enough to detect a significant fraction of the relatively faint submm 
sources. As can be seen for the specific submm galaxies shown in Figs. 14-20, 
when they are available, deep high-resolution radio images are very useful for 
determining accurate positions and even astrophysical properties of submm 
galaxies: see Ivison et al. (2001). 

3.2.3 CO rotation line emission and continuum mm-wave inter jerometry 

The detection of CO line emission from submm-selected galaxies is a crucial 
step in the confirmation of their identification. It has been demonstrated in 
only three cases so far (see Figs. 14, 18 and 19), using the OVRO Millimeter 
Array (MMA; Frayer et al., 1998, 1999), in one case in combination with the 
BIMA array (Ivison et al., 2001), and the IRAM PdBI (Kneib et al., 2002). 
These observations are very time-consuming, typically requiring tens of hours 
of observing time. The ~ 1 GHz instantaneous bandwidth of existing line- 
detection systems also means that a redshift accurate to at least 0.5% must be 
known before attempting a CO detection. In other cases, continuum emission is 
detected using the interferometers, confirming the reality of the initial submm 
detection and providing a better position (Downes et al., 1999a; Bertoldi et al., 
2000; Frayer et al., 2000; Gear et al., 2000; Lutz et al., 2001; Dannerbauer et al., 
2002), but no absolute confirmation of a correct optical/near-IR identification 
or a crucial redshift. 

The ALMA interferometer array will have the collecting area and bandwidth 
to make rapid searches for CO line emission in the direction of known submm 
continuum sources from about 2010. Specialized wide-band mm-wave spectro- 
graphs to search for multiple high- redshift CO lines separated by 115 GHz/ (1 + 
z) that are currently under development (Glenn, 2001). Wide-band cm-wave 
receivers at the 100-m clear- aperture Green Bank Telescope (GBT) could de- 
tect highly redshifted 115-GHz CO(l — > 0) line emission. 

3.2.4 X-ray observations 

Based on synthesis models of the X-ray background radiation intensity (Fabian 
and Barcons, 1992, Hasinger et al., 1996), Almaini et al. (1999) and Gunn and 
Shanks (2002) suggested that 10-20% of the submm galaxy population could 
be associated with the hard X-ray sources that contribute this background. 
Observations of fields with common deep Chandra and SCUBA images were 
discussed in Section 2.8. 

The small degree of observed overlap between the submm and X-ray sources, 
implies that if a significant fraction of submm galaxies are powered by ac- 
cretion in AGN, then the accretion must occur behind an extremely thick 
absorbing column of gas, and less than 1% of the X-ray emission from the 
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AGN can be scattered into the line of sight (Fabian et al., 2000; Barger et 
al., 2001; Almaini et al., 2002). In order to avoid detection using SCUBA, 
high-redshift hard X-ray Chandra sources must either contain a very small 
amount of gas and dust, and thus have only a small fraction of their energy 
reprocessed into the far-IR waveband, which seems unlikely given their hard 
spectra; or they must contain dust at temperatures much higher than appears 
to be typical for submm-selected galaxies. The detection of Chandra sources 
in Abell2390 using ISO at 15 /zm, but not using SCUBA at 850 /im, argues 
in favor of at least some Chandra sources having very hot dust temperatures 
(Wilman et al., 2000). Comparison of larger, deep ISO 15- /im images with 
Chandra images shows that most of the faint, red AGN detected by Chandra 
are detected in the mid-IR (Franceschini et al., 2002). 

This should be readily confirmed using wide-field sensitive mid-IR observations 
of Chandra fields using SIRTF, images which will also yield well-determined 
SEDs for the dust emission from the detected galaxies. 

3.2.5 Mid- and far-IR observations 

Distant submm galaxies are too faint at far- and mid-IR wavelengths to have 
been detected in the all-sky IRAS survey. However, the first submm-selected 
galaxies were detected while the next-generation ISO space observatory was 
still operating, and there were both late-time ISO observations of submm 
fields, and some serendipitous overlap of fields. In general, the small aperture 
and small-format detector arrays of ISO still led to relatively little overlap 
between SCUBA and ISO galaxies, for example in the HDF-N (Hughes et al., 
1998; Elbaz et al., 1998) and Abell2390 (Fabian et al, 2000; Wilman et al., 
2000). The brightest SCUBA galaxies in the cluster Abell370 (Ivison et al., 
1998a) have fluxes measured by ISO at 15 /im (Metcalfe, 2001), providing valu- 
able, and otherwise difficult to obtain, constraints on their short-wavelength 
SEDs (Fig. 2), and thus their dust temperatures. The much larger detector ar- 
rays aboard SIRTF should allow many more constraints to be imposed on the 
mid-IR SEDs of submm galaxies. Existing deep submm fields are included for 
imaging within the SIRTF guaranteed time programs, and individual submm 
galaxies have been targeted for mid-IR spectroscopy.^ 

3.3 A gallery of follow-up results 

In this section we show some of the submm-selected galaxies with the most 
complete and comprehensive follow-up information, including all three with 
confirmed redshifts (Figs. 14, 18 and 19). To reveal more of the diversity of 

13 Details of SIRTF observing programs can be found at sirtf.caltech.edu. 
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counterparts to the SCUBA galaxies, we also show a relatively strong radio 
source with a faint red iC-band counterpart (Fig. 17), two K ~ 19.5 galaxies, 
one a formal ERO and the other with very red colors (Smail et al., 1999; 
Figs. 16 and 20), which are likely to be the correct counterpart on the grounds 
of the relatively low surface density of EROs; and a mm-continuum source 
located by the OVRO MMA at the position of a SCUBA-selected galaxy, 
with a very faint K-band counterpart observed using the NIRC instrument on 
the Keck telescope (Frayer et al., 2000; Fig. 15). Similar enigmatic faint red 
galaxies have been reported as counterparts to well-located submm-selected 
galaxies from other submm surveys by Gear et al. (2000) and Lutz et al. (2001), 
while Dannerbauer et al. (2002) do not find counterparts to three well-located 
MAMBO galaxies to a 3-cr limit of K s = 21.9. 

The galaxies shown in Figs. 14-20 are certainly an unrepresentative sample 
of submm-selected galaxies, missing galaxies that are either intrinsically very 
faint at other wavelengths or lie at the highest redshifts. It is to be hoped 
that within the next few years, deep follow-up observations, especially near- 
IR ground-based observations^] and mid-IR observations using SIRTF will 
reveal the nature of the majority of submm-selected galaxies-p 5 ] 

3.4 Clustering properties 

The deep submm-wave surveys made to date provide relatively little infor- 
mation about the spatial distribution of the detected galaxies. There is some 
indication from the UK 8-mJy SCUBA survey (Almaini et al., 2002; Fox et al., 
2002; Ivison et al., 2002; Scott et al., 2002) and from the widest-field MAMBO 
surveys (Carilli et al., 2001) that the clustering strength of SCUBA-selected 
galaxies is greater than that of faint optically-selected galaxies, yet less than 
that of K ~ 20 ERO samples (Daddi et al, 2000). Webb et al. (2002b) point 
out that the angular clustering signal expected in submm surveys is likely to 
be suppressed by smearing in redshift, as the submm galaxies should have a 
wide range of redshifts, and so the spatial correlation function might in fact be 
stronger than that of the EROs. As most K ~ 20 EROs are evolved elliptical 
galaxies at z ~ 1, which are likely to be amongst the first galaxies to form 
in the most overdense regions of the Universe, their strong clustering is easily 
explained. However, a definitive result on the clustering of submm galaxies 
awaits a much larger sample of galaxies than considered by Scott et al. (2002) 
and Webb et al. (2002b). Characteristic brightness fluctuations on the angular 
scales expected from faint unresolved dusty galaxies have been found in both 

1 Near-IR imaging to 2a K ~ 24 of all of the SCUBA Lens Survey submm galaxies 
is underway using the NIRC camera at the Keck telescopes: see Fig. 15. 
15 More information about this sample can be found in the catalog paper of the 
SCUBA Lens Survey (Smail et al., 2002 and references therein). 
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Fig. 14. Multi-waveband images of SMM J02399— 0136 (23mJy; Ivison et al., 1998; 
Frayer et al., 1998). The format of this figure is the template for those that follow. 
Note that these multiwaveband figures are presented in order of reducing 850-/xm 
flux density, without correcting for gravitational lensing amplification. The leftmost 
panel shows black contours of 850-/xm emission superimposed on a grayscale i-band 
image. The second panel shows black contours of faint 1.4-GHz radio emission su- 
perimposed on a If -band image. These two left-hand images are both 30arcsec 
on a side. The third panel shows a 10-arcsec zoom of the if -band image (from 
UKIRT unless otherwise stated; Smail et al., 2002). The rightmost panel shows a 
-B-band CFHT image in this figure; in the figures that follow this panel shows an 
HST image. Here and in the figures that follow, white contours are added to show 
contrast in saturated regions of the grayscale. North is up and East is to the left. 
SMM J02399— 0136 is a merging galaxy with a confirmed optical/radio counterpart, 
and a CO redshift z = 2.808: see Vernet and Cimatti (2001) for a new high-quality 
spectrum showing Lyman-a emission from this galaxy extended over 12 arcsec. 




Fig. 15. Images of SMM J00266+1708 (18.6 mJy; Frayer et al., 2000). The left-hand 
if -band image is from UKIRT; the right-hand if -band image is from Keck-NIRC. 
The if -band detection is located at the position of the very red galaxy M12 in a 
1.1-mm continuum image obtained using the OVRO MMA. 




Fig. 16. Images of SMM J0942+4658 (17.2 mJy), an ERO counterpart (Smail et al., 
1999). Faint radio emission and extended, rather bright if -band emission make this 
a good candidate for the source of the submm emission. HI is a low-redshift spiral 
galaxy in the foreground of Abell851. 
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Fig. 17. Multi-waveband images of SMM J14009+0252 (14.5 mJy), the bright ra- 
dio-detected submm galaxy Abell 1835 (Fig. 1; Ivison et al., 2000). Two faint near-IR 
counterparts can be seen in the if-band image. Of these, J5 is extremely red, has no 
counterpart in the HST-F702W image, and is aligned accurately with the centroid 
of the radio emission. 




Fig. 18. Multi-waveband images of SMM J14011+0252 (12.3 mJy; Ivison et al., 2000, 
2001). This complex merging system has a confirmed optical/radio counterpart, and 
a CO redshift z = 2.565 (Frayer et al., 1999). High-resolution CO and radio images 
are presented in Ivison et al. (2001). Note that the Northern extension of Jl is 
extremely red, and is close to the centroid of the radio emission. J2 is blue, while 
Jl is red. The complexity of this system is a caution against simple treatment of 
extinction as a uniform screen in submm galaxies: for a detailed discussion see 
Goldader et al. (2002) and references therein. 




Fig. 19. Multi-waveband images of SMM J02399-0134 (11.0 mJy; Kneib et al., 
2002). This ring galaxy has a confirmed optical/radio counterpart, and a CO redshift 
z = 1.06. Its low redshift accounts for its very bright i^-band image and mid-IR 
ISO detection at 15 /um. The other galaxy in the K-band image is a member of 
Abell 370. 



the 850-/xm SCUBA image of the HDF-N (Peacock et al., 2000) and in deep, 
confused 175-//m ISO images (Lagache and Puget, 2000; Kiss et al., 2001). 

Haiman and Knox (2000) have discussed the details of measurements of the 
correlation function of unresolved submm galaxies on arcminute angular scales 
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Fig. 20. Images of SMM J04431+0210 (7.2 mJy), a very red counterpart (Smail 
et al., 1999). A tentative Ha redshift of z = 2.5 is determined from a near-IR 
Keck-NIRSPEC observation. Unlike SMM J09429+4658 this galaxy has no radio 



emission. 



in the context of CMB experiments, finding that the correlated signal can carry 
important information about the nature and evolution of the submm galaxy 
population. A simpler investigation by Scott and White (1999) drew similar 
conclusions, while there is further discussion by Magliocchetti et al. (2001). 



4 Submm galaxy luminosity functions and their relationship with 
other populations 



Submm-selected galaxies are an important component of the Universe, but 
are typically very faint in other wavebands, and so difficult to study. This 
immediately implies that the submm population does not overlap significantly 
with other types of high-redshift galaxies, although it may be possible to infer 
their properties if detailed information about these other classes is available 
(Adelberger and Steidel, 2000). The likely lack of overlap is reinforced by 
the relatively low surface density of submm galaxies as compared with the 
faintest optically-selected galaxies. Nevertheless, a detailed understanding of 
the process of galaxy formation demands that the relationship of the submm 
galaxies to other populations of high-redshift galaxies is determined. 



4-1 Optically- selected Lyman-break galaxies (LBGs) 



The LBGs (Steidel et al., 1999) are sufficiently numerous to have a well- 
defined luminosity function (Adelberger and Steidel, 2000). The effects of dust 
extinction on the inferred luminosity of a small subset of LBGs have been 
estimated reliably from near-IR observations of Ha emission: corrections by 
factors of 4-7 are indicated (Pettini et al, 1998, 2001; Goldader et al. 2002). At 
present, it is difficult to confirm this degree of extinction directly, as attempts 
to detect LBGs using submm-wave instruments have not so far been successful: 
see Section 2.8. Observations suggest that a typical LBG has a 850-/xm flux 
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density of order 0.1 mJy, well below the confusion limit at the resolution of 
existing submm-wave images. 

Adelberger and Steidel (2000) have discussed the various selection effects as- 
sociated with submm, optical and faint radio selection of high-redshift galaxy 
samples. They assumed that the relation between the slope of the UV SED 
of a galaxy and the fraction of its luminosity emitted in the far-IR waveband 
that is observed for low-redshift IUE starbursts with luminosities less than 
about 10 11 L (Meurer et al., 1999) holds at greater redshifts and luminosities. 
A common, smooth luminosity function can then account for the properties 
of LBGs and submm galaxies. A priori, there must be an underlying multi- 
waveband luminosity function of all high-redshift galaxies from which both 
classes of galaxies are drawn. However, while observations of some submm 
galaxies (a key example being SMM J14011+0252; Ivison et al., 2000a, 2001; 
Fig. 18) seem to support this interpretation at first sight, it is clear that only 
the J2 region of this galaxy would be identified as a LBG, while the submm 
emission is concentrated nearer to Jl. Further discussion can be found in 
Goldader et al. (2002). Because of the apparent diversity of optical-submm 
properties of submm galaxies (Ivison et al., 2000a; Smail et al., 2002), this 
simple transformation is unlikely to hold. Hence, a fraction of submm galaxies 
will probably never be detected in rest-frame UV continuum surveys because 
of their extreme faintness. The confirmed ERO submm galaxies (Smail et al., 
1999; Gear et al., 2000; Lutz et al., 2001) are clear examples of such a popu- 
lation. 



4-2 Extremely Red Objects (EROs) 



The development of large format near-IR detectors has enabled relatively deep, 
wide-field IR surveys, and lead to the discovery of a class of faint EROs (galax- 
ies with colors in the range R — K > 5.5 — 6), supplementing traditional low- 
mass stellar EROs (Lockwood, 1970). At first EROs were found one by one, 
(Hu and Ridgeway, 1994; Graham and Dey, 1996), but statistical samples are 
now detected in K < 20 near-IR surveys (Thompson et al., 1999; Yan et al., 
2000; Daddi et al., 2000; Totani et al., 2001), in parallel to their identification 
in multiwaveband surveys (Smail et al., 1999; Pierre et al., 2001; Smith et al., 
2001; Gear et al., 2000; Lutz et al., 2001). The clustering of relatively bright 
K < 19.2 EROs (Daddi et al., 2000) is observed to be very strong, fueling 
speculation that EROs are associated with the deepest potential wells that 
have the greatest density contrast at any epoch. This has been used as an 
argument in favor of their association with submm galaxies, which could be 
good candidates for massive elliptical galaxies in formation (Eales et al., 1999; 
Lilly et al, 1999; Dunlop, 2001). 
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There are two obvious categories of extragalactic EROs: very evolved galax- 
ies, containing only cool low-mass stars, and strongly reddened galaxies, with 
large amounts of dust absorption, but which potentially have a very blue un- 
derlying SED. Only the second are good candidates for identification with 
submm-luminous galaxies. Detections of faint radio emission associated with 
young supernova remnants in EROs, and determinations of signatures of ongo- 
ing star-formation in their rest-frame UV colors should allow these cases to be 
distinguished. Radio and submm follow-up observations of EROs have tended 
to show that most are passively evolved non-star-forming galaxies without 
detectable radio emission (Mohan et al., 2002), with at most about 10-20% 
being candidates for dust-enshrouded starbursts/AGN. It is important to re- 
member that few EROs selected from wide-field near-IR surveys, which reach 
limits of K ~ 20, are actively star- forming submm galaxies (see the summary 
of results in Smith et al., 2001). 

A small but significant fraction of submm galaxies appear to be associated with 
EROs at bright magnitudes K < 20 (Smail et al., 1999, 2002). Many more 
submm galaxies probably fit the ERO color criterion, but at much fainter 
magnitudes; see for example the K = 22.5 SCUBA galaxy shown in Fig. 15 
(Frayer et al., 2000). It is certainly possible that future fainter ERO samples 
with K > 22 could contain a greater fraction of submm-luminous galaxies and 
fewer passive ellipticals than the K ~ 20 samples. Note also that EROs have 
unfavorable K corrections for detection at high redshifts (Dey et al., 1999; 
Gear et al., 2000): beyond z ~ 2.5 any ERO would be extremely faint at even 
near-IR wavelengths. This could account, in part or in whole, for the extreme 
faintness of counterparts to submm galaxies, if many do have extremely red 
intrinsic colors. 



4-3 Faint radio galaxies 

As discussed above, the faintest radio galaxies should be detectable in the 
submm waveband if the far-IR-radio correlation remains valid at high red- 
shifts. The narrow dispersion of this correlation suggests that submm galaxies 
and faint radio galaxies are perhaps the most likely populations of high-redshift 
galaxies to overlap substantially. Surveys made using SCUBA to search for 
optically-faint, and thus presumably high-redshift, galaxies with radio flux 
densities close to the detection threshold of the deepest radio surveys (Barger 
et al., 2000; Chapman et al., 2001b) have been used to detect many tens of 
high-redshift dusty galaxies much more rapidly than blank-field surveys. The 
selection effects at work when making a radio-detected, optically-faint cut 
from a radio survey are not yet sufficiently well quantified to be sure that 
these catalogs are representative of all submm galaxies. The typical optical 
magnitudes of the radio-selected objects with submm detections are clustered 
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around / ~ 24 and greater. Hence, bright optical counterparts to mJy-level 
submm galaxies are rare (Chapman et al., 2002b). Because relatively accurate 
positions are available from the radio observations, it should be possible to 
determine spectroscopic redshifts for a significant number of these galaxies, 
providing a valuable contribution to our knowledge of the distances to at least 
a subset of submm galaxies. 

4-4 Active galaxies and X-ray sources 

An important category of objects that could be associated with submm galax- 
ies are accreting AGN. As discussed in Sections 2.8 and 3.2.4, the overlap be- 
tween Chandra and SCUBA galaxies does not appear to be very great — at 
about the 10% level. However, the gravitational energy released when forming 
the supermassive black holes in the centers of galactic bulges (Magorrian et 
al., 1998) is likely to be only a few times less than the energy released by 
stellar nucleosynthesis over the lifetime of the stars in the bulge, and so a case 
remains for the existence of both a significant population of Compton-thick 
AGN submm sources with no detectable X-ray emission at energies less than 
10 keV, and for hot dusty AGN detectable in the mid- and far-IR but not the 
submm wavebands (Wilman et al., 2000; Blain and Phillips, 2002). 

4-5 Gamma-ray burst (GRB) host galaxies 

An interesting new development is the idea that if GRBs are likely to be asso- 
ciated with the deaths of massive stars, then the rate of GRBs and the global 
high-mass star- formation rate should be linked (Krumholz et al., 1998). By 
searching for submm emission from the directions of GRBs, it may be pos- 
sible to test whether either submm or UV-bright galaxies are the dominant 
population to host high-mass stars, and what fraction of the submm galax- 
ies are powered by non-GRB-generating AGN (Blain and Natarajan, 2000). 
About 10% of GRBs are expected to be in hosts with 850-/im flux densities 
greater than 5mJy (Ramirez-Ruiz et al., 2002), if submm galaxies dominate 
the cosmic star-formation rate and are not typically powered by AGN. Two 
excellent candidates for submm-loud host galaxies of GRBs are now known 
(Berger et al., 2001; Frail et al., 2002), the first based on deep VLA radio im- 
ages, the second on direct SCUBA and MAMBO mm/submm observations. 
Most GRB hosts appear to be associated with R ~ 25 optical galaxies (Bloom 
et al., 2002), which could also be typical of the submm galaxy population. It is 
difficult to detect GRB host galaxies without hitting the confusion limit using 
SCUBA, but attempts are underway. As a byproduct of surveys for submm 
afterglow emission, Smith et al. (1999, 2002) imposed limits to the submm 



53 



host galaxy emission from the direction of 12 GRBs. An ongoing JCMT pro- 
gram (Nial Tanvir et al.) is searching directly for submm emission from the 
host galaxies of accurately located GRBs (Barnard et al., 2002). Detecting 
and resolving the submm emission from GRB host galaxies should ultimately 
be very simple using ALMA, requiring observations of only a few minutes 
minutes per burst. 

4-6 Prospects for the follow-up observations in the future 

In order to make detailed submm-wave studies of the astrophysics of high- 
redshift galaxies, high-resolution images will be required. Existing mm-wave 
interferometers can provide high-quality images of the brightest submm galax- 
ies (for example Frayer et al., 2000; Gear et al., 2000; Lutz et al., 2001); how- 
ever, in order to detect rapidly typical LBGs, EROs and hard X-ray sources, 
the additional sensitivity and resolution that should be provided by ALMA 
is required. Improved information on the SEDs of dusty galaxies will be ob- 
tained using the SIRTF, SOFIA, ASTRO-F and Herschel space- and air-borne 
observatories: see Tables 1 and 2. 



5 Modeling the evolution of submm galaxies 

As soon as the first submm galaxies were detected in 1997, it was clear that 
they made a significant contribution to the luminosity density in the high- 
redshift Universe, subject to the plausible assumptions that their SEDs were 
similar to those of luminous low-redshift dusty galaxies, and that their red- 
shifts were not typically less than about z = 0.5. These assumptions remain 
plausible and have been confirmed to an acceptable level by subsequent ob- 
servations (Smail et al., 2000, 2002). Despite an initial suggestion that 30% 
of the SCUBA galaxies could be at z < 1 (Eales et al., 1999), it now seems 
that a median redshift of submm-selected galaxies is of order 2-3 (Eales et al., 
2000; Smail et al., 2002). 

We have already discussed that evolution by a factor of about 20 in the value 
of L* is required to account for the properties of the submm source counts and 
background radiation intensity. That is the key result from submm surveys, 
but how can it be explained in models of galaxy evolution? In this section we 
will not describe the modeling process in great detail, but we highlight the key 
features of such analyses, and the most important future tests of our current 
understanding. 

It is important is to be aware that there is still considerable uncertainty in the 
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exact form of evolution required to explain the submm observations. While 
strong luminosity evolution of the dusty galaxy population is required out to 
z ~ 1 and beyond, the detailed form of that evolution is rather loosely con- 
strained by count and background data. Redshift distributions are an essential 
requirement in order to determine the form of evolution accurately. 



5.1 An array of possible treatments 

A variety of approaches have been taken to making predictions for and inter- 
preting the results of submm surveys. This work began after the deep 60-/xm 
counts of IRAS galaxies were derived, and it was clear that strong evolution 
was being observed out to z ~ 0.1 (Hacking and Houck, 1987; Saunders et 
al., 1990; Bertin et al., 1997). Observations of more distant galaxies at longer 
wavelengths could probe the extrapolated form of evolution, and disentangle 
the degenerate signatures of density and luminosity evolution (Franceschini et 
al., 1988; Oliver et al. 1992). 

Before the first deep submm survey observations in 1997, Franceschini et al. 
(1991), Blain and Longair (1993a, b, 1996) and Pearson and Rowan-Robinson 
(1996) made a variety of predictions of what might be detected in submm 
surveys. Guiderdoni et al. (1998) and Toffolatti et al. (1998) did similarly 
as the first observational results become available. Generally, the observed 
surface density of submm galaxies was underpredicted. Prior to the recognition 
that an isotropic signal in the far-IR COBE-FIRAS data was an extragalactic 
background (Puget et al., 1996) and not Zodiacal emission (Mather et al., 
1994), this could be accounted for by the use of this unduly restrictive limit 
on the intensity of the submm background spectrum. Once submm count 
data became available after SCUBA was commissioned, it was possible to 
take either a more empirical or a more theoretical view of the consequences 
for dusty galaxy evolution. 

On the empirical side, forms of evolution of the low-redshift luminosity func- 
tion of dusty galaxies (Saunders et al., 1990; Soifer and Neugebauer, 1991) 
that were required to fit the count and background data could be deter- 
mined (Malkan and Stecker, 1998, 2001; Blain et al., 1999b; Tan et al., 1999; 
Pearson, 2001; Rowan-Robinson, 2001). At redshifts less than about unity, 
this is done by requiring that the predicted counts of low-redshift IRAS and 
moderate-redshift ISO galaxies are in reasonable agreement with observations. 
At greater redshifts, where the form of evolution is constrained by submm 
count and background data, there is significant degeneracy in the models: 
strong evolution could proceed all the way out to a relatively low cutoff red- 
shift (z ~ 2-3), or the strong evolution could terminate at a lower redshift 
z ~ 1, followed by a tail of either non-evolving or declining luminosity den- 
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sity out to greater redshifts (z > 5) (see Fig. 9 of Blain et al., 1999b). This 
degeneracy occurs because the far-IR background radiation (like almost all 
backgrounds) is generated predominantly at z ~ 1, while submm galaxies can 
contribute to the counts equally at almost any redshift 1 < z < 10. It can be 
broken by determining a redshift distribution of submm galaxies, which would 
be very different in the two cases. The form of evolution that is consistent with 
the latest observational constraints and radio-derived redshift information for 
submm-selected galaxies (Smail et al., 2002) is shown by the thick solid and 
dashed lines in Fig. 21. 

Note that the assumptions that underlie these derivations are not yet all ver- 
ified by observations. It is unclear whether all high-redshift dusty galaxies 
detected in submm surveys have similar SEDs. It is possible that the proper- 
ties of the dust grains in galaxies evolve with redshift, leading to a systematic 
modification to the temperature or emissivity index. It is reasonable to ex- 
pect the dust-to-gas ratio in the highest redshift galaxies to be lower than in 
low-redshift galaxies, as less enrichment has taken place. However, note that 
enrichment proceeds very rapidly once intense star formation activity is un- 
derway. Even the very first regions of intense star formation could thus be 
readily visible in the submm, despite the global metallicity being extremely 
low. While it seems unlikely, based on a handful of observations (Fig. 2), it 
is certainly possible that a population of dusty galaxies with a significantly 
different SED is missing from current calculations (Blain and Phillips, 2002). 

A more theoretically-motivated approach, based on making assumptions about 
the astrophysical processes at work in galaxy evolution and then predicting the 
observational consequences, has rightly become popular in recent years. These 
'semi-analytic' models, which were generally developed to explain optical and 
near-IR observations, take a representative set of dark-matter halos that evolve 
and merge over cosmic time, from the results of N-body simulations, and de- 
termine their star-formation histories and appearance using a set of recipes for 
star- formation and feedback (White and Frenk, 1991; Kauffmann and White, 
1993; Cole et al, 1994, 2000; Guiderdoni et al., 1998; Granato et al, 2000, 
2001; Baugh et al., 2001; Benson et al., 2001; Somerville et al., 2001). Unfortu- 
nately, at present there is insufficient information from submm observations to 
justify a model that contains more than a handful of uncertain parameters, and 
so it is difficult to exploit the full machinery of semi-analytic models to explain 
the submm observations. Despite the free parameters available, semi-analytic 
models have had limited success in accounting for the observed population of 
high-redshift submm galaxies, without adding in an extra population of more 
luminous galaxies to the standard prescription (Guiderdoni et al., 1998), or 
breaking away from their traditional reliance on a universal initial mass func- 
tion (IMF). As more information becomes available, then the full capabilities 
of the semi-analytic models can hopefully be applied to address dusty galaxy 
evolution. 
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Fig. 21. The history of energy generation in the Universe, parameterized as a 
star-formation rate per unit comoving volume. The absolute normalization of the 
curves depends on the assumed stellar IMF and the fraction of the dust-enshrouded 
luminosity of galaxies that is generated by AGN. The points show results derived 
from a large number of optical and near-IR studies, for which detailed references can 
be found in Blain et al. (1999b) and Smail et al. (2002). The most important results 
are from Lilly et al. (1996; filled stars) and Steidel et al. (1999; high-redshift diagonal 
crosses). The up-pointing arrow comes from the submm-based estimate of Hughes 
et al. (1998). An important new measurement of the extinction- free low-redshift 
star-formation rate from radio data, that is not plotted, has been obtained by Yun 
et al. (2001): 0.015 ± 0.005 M Q yr _1 . The thick solid and dashed lines represent the 
current best fits to far-IR and submm data in a simple luminosity evolution model 
and a hierarchical model of luminous merging galaxies respectively, as updated to 
reflect additional data and a currently favored non-zero-A cosmology. The thinner 
solid lines show the approximate envelope of 68% uncertainty in the results of the lu- 
minosity evolution model. The thin and thick dotted lines represent the best-fitting 
results obtained in the original derivations (Blain et al., 1999b, c). 

In a blend of these approaches, assuming that the SCUBA galaxies are all 
associated with merging galaxies (Ivison et al., 1998a, 2001; Figs. 14 and 18), 
and yet not being sure of the physical processes by which the mergers generate 
the luminosity we observe, Blain et al. (1999c) used a minimally-parametrized 
semi-analytic model to investigate the change in the properties of merging 
galaxies required to reproduce the submm and far-IR counts and backgrounds: 
see also Jameson (2000) and Longair (2000). The observed background spec- 
trum can only be reproduced for strong evolution of the total luminosity den- 
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sity out to redshift z ~ 1, by a factor of about 20 (see Fig. 21). In addition, the 
lifetime of the luminous phase associated with mergers, and thus the mass-to- 
light ratio also had to be reduced by a large factor at high redshift in order 
to reproduce the observed submm counts. The physical reason for this change 
must be an increased efficiency of star formation during starburst activity or 
of AGN fueling at increasing redshifts; both make sense in light of the greater 
gas densities expected at high redshifts. The model has the advantage of being 
able to reproduce the faint optical counts, if blue LBGs are also associated 
with merging galaxies. The submm galaxies release about four times more en- 
ergy in total than the LBGs, and do so over a period of time during a merger 
that is about ten times shorter. It is likely that the total baryonic mass and 
geometry of the merging galaxies also play important roles in determining the 
details of star-formation activity or AGN fueling during a merger. 

When reviewing the predictions and results of any model, note that it is easy to 
produce a model that can account for the far-IR-submm background radiation 
intensity; more difficult to account for the submm counts; and more difficult 
again to reproduce a plausible redshift distribution. 



5.2 Observational tests of models 

The key observational test of models of submm-wave galaxy formation is the 
redshift distribution, which is know in outline from radio-submm observa- 
tions (Smail et al., 2000). Determining the redshift distribution is a key goal 
of extensive ongoing follow-up observations, but the process has proved to 
be difficult and time-consuming, as documented extensively by Smail et al. 
(2002). The crucial problems are the faintness of the counterparts, combined 
with the relatively poor positional accuracy of the centroids of the submm 
galaxies, which are unresolved due to the coarse spatial resolution of existing 
submm images. 

Detailed measurements of the counts of galaxies at both brighter and fainter 
flux densities than those shown in Fig. 9 would also constrain models. However, 
determining the bright counts requires a large-area survey, which is likely to 
be relatively inefficient (Fig. 23), while determining the faint counts requires 
greater angular resolution than can be provided by the telescopes used to 
make existing surveys, to avoid source confusion. The very bright counts will 
certainly be probed directly towards the end of the decade by the Planck 
Surveyor all-sky survey at a resolution of 5 arcmin, and sooner by large-area 
surveys using forthcoming large-format mm/submm-wave bolometer arrays 
on ground-based telescopes, including BOLOCAM (Glenn et al., 1998) and 
SCUBA-II. Limits on the bright submm-wave counts can be imposed from 
the number of candidate point sources that can be found in large-area submm 
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maps of Galactic fields (Pierce-Price et al., 2001; Barnard et al., 2002). The 
faint counts will ultimately be determined directly using the SMA, CARMA 
and ALMA interferometers. 

The results of deep ISO surveys have been regularly cited as a useful constraint 
on galaxy evolution (Rowan-Robinson et al., 1997; Xu, 2000; Chary and El- 
baz, 2001). This is certainly true out to z ~ 1. However, when estimating a 
total luminosity density from 15-/xm data, it is vital that the correct SED is 
used to extrapolate to longer wavelengths, as it is easy to overestimate the 
amount of luminosity associated with a 15-/zm source by assuming a mid-IR 
SED that is too steep. For example, compare the inferred luminosity density 
results at redshifts z ~ 0.7 quoted by Rowan-Robinson et al. (1997) and Flores 
et al. (1999). The results differ by a factor of 5; Flores et al. (1999) obtain the 
lower result by using radio observations to constrain the total luminosity of 
the galaxies detected at 15 /zm. Extrapolating mid-IR data towards the peak 
of the SED at longer wavelengths is more difficult than extrapolating submm 
observations to fix the position of the peak of the SED that lies at shorter 
wavelengths. This is both because the form of the SED is intrinsically sim- 
pler on the long-wavelength side of the peak, and because the well-determined 
spectrum of the far-IR background radiation can be used to constrain the 
luminosity-averaged dust temperature of the submm galaxies. Mid-IR obser- 
vations with SIRTF after 2002 will provide much more information about the 
SEDs and evolution of dusty galaxies to redshifts z ~ 2. 

In Fig. 22 we show the deep 15-/xm counts predicted by models designed to 
account for the submm data (Blain et al., 1999b, c), updated to the current 
data and cosmology. If the mid-IR SED is chosen appropriately, then the fit 
is quite acceptable. Including PAH emission features or varying the mid-IR 
SED index a has relatively little effect on the result. The same approach 
can be used to estimate the deep cm- wave radio counts. If we assume just 
the form of the radio-far-IR correlation (Condon, 1992), without any fine 
tuning, and a radio SED of the form f v oc z/ -0 6 , then the predicted 8.4-GHz 
counts brighter than lOyuJy, based on the submm-based models are 1.05 and 
0.98arcmin~ 2 respectively; the corresponding power-law indices of the count 
function N(> S) oc S a are a = —1.4 and —1.3 respectively. The results in both 
models match the observed 8.4-GHz 10-/xJy count of 1.01 ±0.14 arcmin -2 with 
a = —1.25 ± 0.2 (Partridge et al., 1997). The reasonable agreement between 
the predictions of the models, which are constrained only by observations 
in the submm and far-IR, and the observed deep mid-IR and radio counts 
confirms that the models are reliable. The source confusion estimates shown 
in Fig. 13, which are based on the same models, should thus be reliable over 
a wide wavelength range from about 10 cm to 10 /iin. 
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Fig. 22. A summary of observed (Elbaz et al., 1999) and predicted (Blain et al., 
1999b, c) differential counts of galaxies in the 15-^m ISO band. The model predic- 
tions assume only a power-law SED in the mid-IR, with f u oc z^ -1 ' 95 : no fine-tuning 
with PAH features in the SED is included. The hierarchical model (Blain et al., 
1999c) provides a better fit to the data, but both agree reasonably well with the 
observations. Both models asymptote to the same count at greater flux densities. 
Note that the relative form of the model counts reflects that seen in Fig. 9, with the 
hierarchical model having the steepest rise. 

5. 3 Modeling the detailed astrophysics of the submm galaxies 

It is not possible to separate the modeling of the evolution of the popula- 
tion of submm galaxies fully from studies of the nature of the galaxies them- 
selves. Their luminosities and masses (see Frayer et al., 1999) demand that 
the submm-luminous phase be short-lived as compared with the age of the 
Universe. The observational information for most of the submm galaxies is 
insufficient to be confident that their nature is understood at present. Of 
the galaxies with reliable counterparts, there are three bright Class-2 galaxies 
(Soucail et al., 1999; Ivison et al., 1998a, 2000a, 2001; Vernet and Cimatti, 
2001) with optical redshifts and CO detections (Frayer et al., 1998, 1999; 
Kneib et al., 2002), and a total of four Class-1 galaxies (Smail et al., 1999; 
Bertoldi et al., 2000; Gear et al., 2000; Lutz et al. 2001), which are known to 
be either very red galaxies or formal EROs. Other submm-selected galaxies 
with accurate positions from radio observations (Smail et al., 2000) or mm- 
wave interferometry (Downes et al., 1999b; Bertoldi et al., 2000; Frayer et al., 
2000; Dannerbauer et al., 2002) remain enigmatic. All that can be said about 
these galaxies is that they all appear to have thermal dust spectra, are all very 
faint at optical wavelengths, and most also appear to be very faint at near-IR 
wavelengths. 
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The Class-2 galaxies are all clearly undergoing mergers or interactions. Much 
less is known about the morphology of the faint, but typically extended Class- 1 
galaxies (Figs. 16, 17 and 20). It is certainly possible that they too are involved 
in interactions, which appear to trigger the dramatic luminosity of almost all 
the low-redshift ULIRGs (Sanders and Mirabel, 1996). Programs of ultradeep 
near-IR imaging on 10-m-class telescopes should soon test this idea. 

Hydrodynamical simulations of gas-rich mergers by Mihos and Hernquist (1996), 
Bekki et al. (1999) and Mihos (2000) show the formation of very dense concen- 
trations of gas, which could be associated with short-lived, very-intense bursts 
of star formation. However, at present it is not possible to simulate a represen- 
tative sample of mergers with the range of geometries likely to be encountered, 
the necessary time resolution, and a sufficiently accurate treatment of the de- 
tailed astrophysics of star-formation to make a reliable connection between the 
limited observations and the underlying galaxy properties. The spatial extent 
of the three bright Class-2 galaxies in the optical waveband appears to be con- 
siderably greater than that of most low-redshift ULIRGs. It is thus difficult 
to be sure that simulations of well-studied low-redshift ULIRGs adequately 
represent the properties of the high-redshift submm galaxies. Note, however, 
that the precise spatial relationship between the optical and submm emission 
in these objects is still unclear (Ivison et al., 2001); the submm emission could 
be more compact than the optical galaxy. 

Larger submm galaxy samples will be available over the next few years, boost- 
ing the likelihood that examples of the full range of submm galaxies will be 
available to be studied in detail. More sensitive observations of the properties 
of the known galaxies will also improve our knowledge of their astrophysics. 
One key question is the relationship of the submm galaxies to the formation of 
elliptical galaxies (Lilly et al., 1999). Whether the bulk of submm galaxies are 
high-redshift low-angular-momentum gas clouds, forming elliptical galaxies in 
a single episode by a 'monolithic collapse' (Eggen et al., 1962), as advocated 
by for example Archibald et al. (2002), or galaxies observed during one of a 
series of repeated mergers of gas-rich, but pre-existing galaxy sub-units, likely 
to take place at relatively lower redshifts, as discussed by Sanders (2001), and 
which might ultimately yield elliptical merger remnants, is an important ques- 
tion that future follow-up observations will address. Existing observations of 
extended and disturbed counterparts to submm galaxies (Ivison et al., 1998a, 
2001; Lutz et al., 2001) tend to favor the second explanation, in which well- 
defined pre-existing stellar systems merge. However, in either scenario it is 
likely that the bulk of the stellar population in the resulting galaxies form 
during the submm- luminous phase. 
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5.4 The global evolution of dust- enshrouded galaxies 



Fig. 21 summarizes the current state of knowledge of the strong evolution 
of the comoving luminosity density contributed by luminous far-IR galaxies, 
whose emission is redshifted into the submm (Blain et al., 1999b, c; Smail et 
al., 2002). The derivation of these results was discussed briefly above, and is 
explained in much more detail in these papers. The results of both the lumi- 
nosity evolution and hierarchical models, which both include strong luminosity 
evolution of the population of low-redshift IRAS galaxies out to high redshifts, 
are fully compatible with the redshift information available for submm galaxy 
samples. This would not be the case if the redshift evolution of the luminosity 
density was markedly different from the forms shown in Fig. 21. For example, 
if the luminosity density of submm galaxies were to match rather than exceed 
the value denoted by the datapoints at z ~ 1 in Fig. 21, and remain at the 
same high level out to z ~ 10, then the submm-wave spectrum of the cosmic 
background radiation (Fig. 11) would tend to be too flat, the intensity of the 
far-IR background radiation would fall short of the observed level, and the 
predicted redshift distribution of submm galaxies would be biased strongly 
to the highest redshifts, which at present seems not to be the case (Smail et 
al., 2002). The determination of complete redshift distributions for existing 
samples of submm galaxies, and of more redshifts for individual very lumi- 
nous galaxies drawn from large, future submm-selected galaxy catalogs, will 
ultimately allow the evolution of distant dusty galaxies to be traced in detail. 



6 Gravitational lensing in the submm waveband 

The first submm-wave surveys for distant galaxies exploited both the weak to 
moderate gravitational lensing magnification, by about a factor of 2-3, expe- 
rienced throughout the inner few square arcminutes of rich foreground cluster 
of galaxies at a moderate redshift in the range z ~ 0.2-0.4 (Blain, 1997) and 
the greater magnification produced along critical lines for much smaller ar- 
eas of the background sky. A 5-arcmin 2 SCUBA field centered on a moderate 
redshift cluster includes both these regions, enhancing the flux density from 
all high- redshift background galaxies (Smail et al., 1997). More, and in some 
cases deeper, SCUBA images of clusters have been taken (Smail et al., 2002; 
Chapman et al., 2002a; Cowie et al., 2002), especially in Abell2218, where a 
multiply-imaged source has been detected (van der Werf and Kraiberg Knud- 
sen, 2001). Whether the magnification acts to increase the surface density of 
background galaxies on the sky depends on the form of their counts. Lensing by 
both galaxies and clusters could have significant applications in future submm 
surveys, especially those sampling the steep counts of bright submm galaxies 
in wide fields, including the all-sky survey from Planck Surveyor (Blain, 1998), 
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and surveys using BLAST, Herschel-SPIRE and SCUBA-II (Tables 1 and 2) 
covering many tens of square degrees. 

6.1 Magnification bias 

Because surface brightness is conserved by all gravitational lenses, the net 
effect of magnifying a population of background galaxies depends on the slope 
of the counts diV(> S)/dS, where N(> S) is the number of galaxies per 
unit area on the sky brighter than S (Schneider et al., 1992). Subject to a 
magnification //, the count becomes [l//x 2 ]diV[> (S/n)]/dS. For a power-law 
count with dN(> S)/dS oc S a , a value of a < — 2 corresponds to an increase 
in surface density if /x > 1. This threshold value corresponds to a slope of 
— 1 for the integral counts N(> S) shown in Fig. 9. Note that for a uniform 
non-evolving population of galaxies a = —2.5. 

As shown in Fig. 9, submm-wave counts are expected to be steep, and to 
change slope sharply at mJy flux density levels, as compared with deep optical 
or radio counts. The significant changes in the count slope are particularly 
unusual, and not found in any other waveband. As a result, the magnification 
bias can be large, increasing the number of detectable bright galaxies (Blain, 

1996, 1997), and providing a way to investigate very faint counts by comparing 
lensed and unlensed fields; for example in the innermost regions of clusters of 
galaxies (Blain, 2002). 

That a significant magnification bias can be exploited using relatively weak 
lensing by clusters of galaxies can be seen by comparing the number of 10- 
20-mJy 850-/im galaxies detected in the SCUBA Lens Survey (Smail et al., 

1997, 2002), and in the larger field of the unlensed 8-mJy survey (Scott et al., 
2002). The ratio is about 3:1, showing a clear positive magnification bias, and 
indicating that if the 850-/xm counts at flux densities greater than 10 mJy are 
represented by a power-law, then the index a < —2. 

6.2 Conditions for exploiting submm lensing by galaxies 

The key advantage of observing background galaxies that are gravitationally 
lensed by foreground mass concentrations in the submm waveband is that the 
K correction (Fig. 4) acts to brighten the distant background lensed galaxy as 
compared with the lens. This is already very familiar from surveys of lensed 
radio AGN (Rusin, 2001), and is illustrated clearly in Fig. 1, in which only the 
central cD galaxy in the lensing cluster shows any significant submm emission. 

In SCUBA cluster lens surveys, both the image separations, and the extent 
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of the high-magnification regions are of order 1 arcmin, a scale which is well 
matched both to the 15-arcsec resolution of the JCMT and to the 2. 5- arcmin 
field of view of SCUBA. The magnification ensures that a significantly greater 
fraction of the submm-wave background radiation intensity is thus resolved 
into detectable galaxies in surveys in the fields of gravitational lensing clusters 
than in even the deepest blank- field surveys (Blain et al., 1999a). However, 
for background sources lensed by galaxies rather than clusters, the relevant 
image separations and the extent of the high-magnification region are only of 
order 1 arcsec, and so cannot be resolved using any single-antenna telescope. 
High-resolution submm observations are required to disentangle lensed and 
unlensed galaxies; this capability will be provided by ALMA (Blain, 2002), 
while pilot studies of should be possible using the CARMA, SMA and IRAM 
PdBI interferometers. The most luminous lensed sources can already be re- 
solved into multiple images using the IRAM mm-wave interferometer (Alloin 
et al., 1997). 

The only caveat for exploiting galaxy-scale lensing is that the source size must 
be small as compared with the area of sky behind the lens that is strongly 
magnified. The intense far-IR and submm emission from low-redshift ULIRGs 
is typically very compact (several hundred pc across; Downes and Solomon, 
1998), and would easily meet this condition; however, there are indications that 
the dust emission from at least some luminous high-redshift submm galaxies 
could extend over scales greater than 10-kpc (Papadopoulos et al., 2001; Chap- 
man et al., 2001a; Lutz et al., 2001; Isaak et al., 2002; Ivison et al., 2001). The 
whole area of sky covered by these galaxies would not then be lensed efficiently 
by an intervening galaxy, although bright knots of emission within them could 
still be magnified by large factors. This concern about lensing efficiency and 
the angular size of distant submm galaxies does not apply to lensing by much 
larger clusters of galaxies, which will always be effective. 



6. 3 Prospects for the lensing studies in the future 

Larger area surveys for brighter samples of luminous dusty galaxies using the 
array of forthcoming ground-based, air- and space-borne instruments, includ- 
ing BOLOCAM, BLAST, SOFIA, SCUBA-II, SI RTF, Herschel and Planck 
Surveyor (see Table 2), should be subject to an enhanced magnification bias. 
High-resolution follow-up observations using ALMA should then yield a large 
sample of strongly magnified high-redshift lensed systems to complement the 
systematically selected CLASS sample of gravitational lensed AGN identified 
at radio wavelengths (Rusin, 2001). These surveys will not be subject to any 
extinction bias due to absorption by dust in the lensing galaxies, and should 
yield a very complete and reliable catalog of up to several thousand lenses 
(Blain, 1998). 



64 



7 Future developments in submm cosmology 

During the last 4 years, the first steps have been taken towards investigating 
the Universe using direct submm-wave surveys. The technologies of the class of 
detectors that made these initial surveys possible are still developing rapidly. 
Many instrumentation projects are underway, which will allow us to increase 
the sizes of samples of distant submm galaxies, and to study known examples 
in more detail; some of their key features are outlined in Tables 1-3. 



7.1 New technologies for instrumentation 

A key technology under development is for bolometers with superconduct- 
ing temperature-sensitive elements, including transition-edge sensors (TESs). 
These are much more stable than the semiconducting thermistors used in ex- 
isting systems, and so can be read out using multiplexed, and therefore much 
simpler, cold electronics. Another advantage of TES devices is that they re- 
quire no bias current, and so need fewer heat-conducting, difficult-to-assemble 
connections to each device. The prototype Fabry-Perot spectroscopic device 
FIBRE, which uses TES bolometers (Benford et al., 2001), was tested suc- 
cessfully at the CSO in May 2001. 

TES devices offer the prospect of increasing the size of the arrays of detector 
elements in mm/submm-wave cameras from of order 100 to of order 10 4-5 , 
providing much larger fields of view. Filled-array detector devices using con- 
ventional semiconducting bolometers are being demonstrated in the SHARC-II 
and HAWC cameras for the CSO and SOFIA, while the SCUBA-II camera 
has a goal of at least a 8 x 8 arcmin field of view — about 25 times greater than 
the field of view of SCUBA — is under development in Edinburgh and is ex- 
pected to integrate large arrays and superconducting bolometers. SCUBA-II 
will supplement its much wider field of view with an enhanced point source 
sensitivity: the same galaxies should be detectable about 8 and 4 times faster 
using SCUBA-II as compared with SCUBA at wavelengths of 850 and 450 /zm 
respectively: see Table 1. A 10-m telescope operating at 850-//m with a 10 5 - 
element detector would have a square field of view about 1 deg on a side, at a 
Nyquist-sampled plate scale: much larger than the 5-arcmin 2 fields of view of 
SCUBA and MAMBO. With such large fields of view, it is not unreasonable 
to survey most of the sky down to the confusion limit of a 10-30-m telescope 
in an observing campaign lasting for several years. 

By combining large numbers of bolometer detectors with dispersive mm-wave 
optics (Glenn, 2001), it should be possible to obtain low- resolution mm-wave 
spectra of galaxies over a very wide band, perhaps 100 GHz, to search for CO 
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and atomic fine-structure line emission from high-redshift submm galaxies 
detected in continuum surveys, while simultaneously carrying out unbiased 
surveys for line-emitting galaxies within the field of view (Blain et al., 2000b). 
Development of several such systems is underway. 

Phase- sensitive heterodyne submm detectors are already very efficient; how- 
ever, only small numbers of these detectors can currently be fabricated into 
an array. Their strength is in very high resolution submm-wave spectroscopy, 
and as sensitive coherent detectors in existing mm-wave interferometers. They 
will be fitted to the forthcoming SMA and CARMA, and will be be exploited 
to the full with the large collecting area of the ALMA array. 

Sensitive arrays of mid- and far-IR detectors should soon be flying, both in 
space aboard SIRTF and ASTRO-F, and in the upper atmosphere, on balloons 
such as BLAST, and aboard SOFIA. Limits to the continuum flux densities 
of the most luminous high-redshift galaxies derived using these facilities, mea- 
sured close to the peak of their SED (see Fig. 2), will provide valuable infor- 
mation about their properties. Spectrographs, both aboard these facilities and 
on ground-based telescopes, will provide detections of and sensitive limits to 
the line radiation from the same objects, providing redshift information and 
astrophysical diagnostics. 



7.2 New telescopes 

At present, a wide range of submm-wave telescopes are available. Single- 
antenna telescopes include the 10.4-m CSO on Mauna Kea, the 10-m Heinrich 
Hertz Telescope (HHT) on Mount Graham, the 15-m JCMT on Mauna Kea, 
the 15-m SEST at La Silla in Chile, the 30-m IRAM telescope in Spain, and 
the 45-m antenna at Nobeyama in Japan. The Large Millimeter Telescope 
(LMT), a 50-m mm-wave telescope is under construction on a 5000-m peak 
near Puebla in Mexico, and it is hoped that the 100-m Green Bank Telescope 
(GBT) in West Virginia can operate at 90 GHz/3 mm during the winter. New 
single-antenna telescopes with large survey cameras have been proposed for 
the excellent submm observing sites at the South Pole and the ALMA site in 
Chile. The Planck Surveyor CMB imaging mission will generate an all-sky map 
in the submm at a resolution of 5arcmin, and the 3.5-m Herschel space tele- 
scope will carry out pointed submm imaging and spectroscopic observations 
of known galaxies, and carry out deep confusion-limited cosmological surveys 
over fields several hundred square degrees in size. Cameras exploiting the 2.5-m 
telescope aboard SOFIA and BLAST and other dedicated ultra-long-duration 
balloon instruments will allow far-IR and submm-wave observations from the 
upper atmosphere. 
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Existing mm-wave interferometers include the 6 x 15-m IRAM PdBI, the 6 x 
10.4-m OVRO MMA, the 10 x 6-m BIMA array at Hat Creek in California 
and the 6 x 10-m Nobeyama Millimeter Array. The 8 x 6-m SMA is under 
construction on Mauna Kea, the first imaging submm-wave interferometer, 
while it is planned to combine 9 of the BIMA antennas with the OVRO MMA 
at a high site in the Inyo Mountains east of Owens Valley in California to form 
CARMA. The international 64 x 12-m ALMA submm interferometer array 
in Chile will provide a tremendous increase in the capability of submm-wave 
spectral line and continuum imaging, providing 10- to 30-/xarcsec resolution, 
and detailed images of even the most distant galaxies. The most luminous 
submm galaxies so far discovered, with 850-/im flux densities of about 25 mJy, 
could be detected at a 10-er significance by ALMA in about a second. Its 
excellent sensitivity and wide 8-GHz instantaneous bandwidth will allow a 
significant fraction of the galaxies detected in deep surveys to be detected 
simultaneously in the continuum and CO rotation lines, providing direct and 
exact redshifts. As the redshifted ladder of CO lines are separated by 115/(1 + 
z) GHz, about 25% of galaxies at z ~ 2.5 will have a CO line lying within the 
8-GHz-wide ALMA band (Blain et al., 2000b). 

7. 3 Future capabilities and progress 

The enhanced capabilities of this array of new facilities is illustrated in Fig. 23. 
The rate at which galaxies can be detected is likely to grow dramatically from 
a few per day at present to many hundreds per hour. Note that the various 
instruments operate at different wavelengths, and so each is most sensitive to 
galaxies at different redshifts and with different luminosities. However, sample 
sizes are certain to increase dramatically, especially when the 10 4 -10 5 galaxies 
that will be detected in the Planck Surveyor all-sky survey are taken into 
account. 

Multiwavelength follow-up observations of all these new submm galaxies are 
likely to remain a time-consuming challenge. However, the likely availability of 
30-m-aperture ground-based optical/near-IR telescopes in the next decades, 
and the extremely deep imaging capability of NGST, should help us to study 
a complete sample of submm galaxies down to luminosities that are only a 
fraction of L*. 



8 Summary: key questions and targets for the future 

The first generation of extragalactic submm-wave surveys have provided an 
important complement to more traditional optical and radio searches for dis- 
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Table 1 

Wavelengths A, sensitivities (as noise equivalent flux density — NEFD), fields of 
view (FOV), and confusion limits due to galaxies and the ISM (in brackets) for 
existing and future ground-based and airborne instruments. 
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An estimate of the speed of a survey down to a chosen depth can be obtained 
by multiplying the FOV by the the square of the NEFD value. The approximate 
extragalactic confusion noise values are the flux density at which there is one brighter 
source per beam (Blain et al., 1998; Fig. 13). This corresponds approximately to 
the width of peak in the non-Gaussian confusion noise distribution (see Fig. 12). 
The expected ISM confusion noise (in brackets) is calculated for a 100-^m surface 
brightness i?o=lMJysr _1 (Helou and Beichman, 1990), and scales as -Bq 5 - Other 
instruments under development, which have not published detailed performance 
estimates include the 350-/xm SHARC-II camera for the CSO. The FOV and NEFD 
values are chosen to provide the correct results for making a fully-sampled image of 
the sky, not measuring the flux from a single galaxy. Updated from Table 1 in Blain 
(1999b). Relevant references are listed in Table 3. TBD: to be decided. 
1 Note that the FOV of MAMBO is expanded by a factor of 3 for the winter of 
2001/2002, with a 117-bolometer detector array. 2 A similar device SIMBA is being 
commissioned at SEST. 
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Fig. 23. The detection rates expected in a variety of forthcoming mm, submm and 
far-IR surveys. The names and wavelengths in microns of the relevant instruments 
are listed: see Tables 1 and 2. The higher- and lower-peaking SCUBA-II curves 
correspond to wavelengths of 450 and 850 fim respectively. References to instrument 
performance for these calculations (Blain and Longair, 1996) can be found in Table 3. 
The surface density of galaxies assumed follows the models of Blain et al. (1999b): 
see Figs. 9 and 10. Curves stop on the right if the surface density is expected to fall 
below a single galaxy on the whole sky. Curves stop at the left when a relatively 
optimistic definition of the 5a confusion noise level for detection is reached (see 
Fig. 13). 

tant galaxies, and discovered a cosmologically significant population of very- 
luminous, high-redshift dusty galaxies. 

We have found that is very hard to study a complete sample of submm galaxies 
at other wavelengths (Smail et al., 2002). The similar experience of other 
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Table 2 

The equivalent to Table 1 for space-borne instrumentation. 
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SPIRIT 

Note that the values listed for Planck Surveyor apply to an all-sky survey. Another 
instrument under development, which has not published detailed performance esti- 
mates is the 50-200-jum sky survey from the Japanese ASTRO-F/IRIS satellite sky 
survey. 

groups involved in deep mm/submm surveys (Barger et al., 1999a; Eales et al., 
2000; Carilli et al, 2001; Scott et al, 2002; Webb et al., 2002b) is reflected in 
the relatively few papers describing the individual multi-waveband properties 
of the almost 200 galaxies detected. The most sensitive follow-up observations 
are required in the near-IR, radio and optical wavebands to identify and study 
them (Frayer et al., 2000; Ivison et al., 2001), that is very faint detection 
thresholds of order 10 /iJy at 1.4-GHz, K ~ 23 and B 3> 26 respectively. 
Much more time has been devoted to multi-waveband follow-up observations 
than was spent on the initial submm detections. Typical examples of the 
submm population can be detected in imaging-mode SCUBA observations in 
about 10 h of integration. However, at least 2h of near-IR observations at 
the 10-m Keck telescope and about 24 h of integration at the VLA are then 
required in order to find likely counterparts to typical submm galaxies. The 
advantage of the VLA radio observations over those at optical and near-IR 
wavelengths is the very large field of view, which allows many galaxies to be 
detected simultaneously. The very brightest optical counterparts to submm 
galaxies can be identified spectroscopically in about 7h of integration using 
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Table 3 

References to instruments listed in Tables 1 and 2. 



Name 



Information 



SCUBA 

MAMBO 

SCUBA-II 

SOFIA 
BOLOCAM 

BLAST 
SMA 
ALMA 
SKA 

Herschel-SPIRE 
Herschel-PACS 
SIRTF-MIPS & IRAC 
Planck Surveyor 
SPECS I SPIRIT 



Holland et al. (1999) 
Kreysa et al. (1998) 

www . j ach . hawaii . edu/ J ACpublic / J C MT/ 
Continuum_observing/SCUBA-2/home.html 

Davidson et al. (1999); sofia.arc.nasa.gov 

Glenn et al. (1998); www-lmt.phast.umass.edu/ins/ 
continuum/bolocam.html 

Devlin (2001); www.hep.upenn.edu/blast 

Ho (2000); sma2.harvard.edu 

Wootten (2001); www.alma.nrao.edu 

www . nfr a. nl/ skai 

www.ssd.rl.ac.uk/spire 

pacs . ster . kuleuven . ac . be 

sirtf.caltech.edu 

astro.estec.esa.nl/Planck 

Mather et al. (1998); space.gsfc.nasa.gov/astro/specs 



4-m class telescopes (Ivison et al., 1998) and higher- quality spectra can be 
obtained in a comparable time using 8-m class telescopes (see the results 
of a 5-h integration using the UVES spectrograph at the European Southern 
Observatory (ESO) VLT by Vernet and Cimatti, 2001). In all cases, identifying 
a plausible counterpart, where possible, is only a first step; finding a redshift for 
these typically faint, red galaxies is much more challenging. In this context, the 
unusual sensitivity of submm surveys to the most distant galaxies is almost a 
drawback, making it very hard to detect a complete sample of submm galaxies 
at other wavelengths. 

Key questions for understanding submm galaxies in the future include: 

What are the properties of typical submm galaxies in other wavebands, and 
what is their relationship to other high-redshift galaxy samples? The submm- 
selected galaxies appear to a diverse mixture of types, including bright merging 
systems (Ivison et al., 1998a, 2000a, 2001), optical QSOs (Kraiberg Knudsen 
et al, 2001), EROs with K < 20 (Smail et al, 1999, 2002; Gear et al., 2000; 
Lutz et al., 2001), and much fainter IR-detected galaxies (Frayer et al., 2000), 
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which may also turn out to have very red colors. It seems that the overlap 
between the 850-/zm submm galaxy population and both the LBGs and faint 
Chandra X-ray sources is small. Note that some of this apparent diversity is 
sure to be due to the very wide redshift distribution of the submm galaxies. 

What is the redshift distribution of the submm galaxies? Models of the evo- 
lution of submm galaxies that do not grossly violate basic observational con- 
straints on the source counts and cosmic background radiation are easy to 
generate. However, it is vital to predict a plausible redshift distribution, with 
only a small fraction at redshifts less than unity, and a probable median red- 
shift of at least 2-3. It is easy to generate a redshift distribution that is biased 
too high. The observational determination of a redshift distribution for a well- 
defined sample of submm galaxies remains a crucial goal. This will be easy 
with ALMA. In the meantime, concerted and time-consuming campaigns of 
optical and near-IR spectroscopy will pay off gradually, while observations of 
cm-wave megamasers and the development of wide-band mm- and cm-wave 
spectrometers may offer alternative routes. The forthcoming (sub)mm inter- 
ferometers CARMA and SMA, and developments of the IRAM PdBI will also 
provide accurate positions and some CO redshifts for submm galaxies. 

What are the details of the astrophysics responsible for the luminosity of the 
submm galaxies? This is very important, as the submm galaxies appear to be 
signposts to some of the most luminous and violent phases of galaxy evolution, 
and could be associated with the formation of the bulk of galactic bulges, 
elliptical galaxies and supermassive black holes (Lilly et al., 1999; Dunlop, 
2001). Whether these galaxies are formed in a single event, or as a series of 
lesser bursts, is a key question for our understanding of the process of galaxy 
formation and evolution. Detailed comparisons of the luminosity derived from 
dust continuum emission, the dynamical mass inferred from molecular line 
profiles, the evolved stellar mass inferred from near-IR observations, and the 
spatial extent of the activity from various high-resolution observations will all 
be important for disentangling the complex astrophysics of these systems. 

When can submm instruments be used to resolve and study high-redshift 
galaxies in detail? This is already practical given enough observing time at 
the OVRO MMA and the IRAM PdBI. The CARMA and SMA interferom- 
eters will soon have important roles to play in these studies. In about 10 
years, ALMA will provide the first real chance to detect and study galax- 
ies rapidly and in great detail using submm observations alone. Luminosities, 
redshifts, dynamical masses and metallicities could all be determined with- 
out needing to resort to radio, optical and near-IR observations as a matter 
of course. However, because ALMA has a relatively small field of view, the 
most efficient survey strategy may be to detect large numbers of galaxies using 
wide-field mm/submm cameras like BOLOCAM, SCUBA-II and their succes- 
sors on single-antenna 10-50-m aperture survey telescopes, and the Herschel 
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and Planck Surveyor space missions, and then use ALMA to provide detailed 
images and spectra of all the detected galaxies. 

What is the fundamental limit to making submm observations of distant galax- 
ies? Submm observations rely on the presence of metals, in the form of molec- 
ular gas or dust grains in order to detect galaxies. While submm radiation 
is able to travel unattenuated across the Universe from prior to the epoch of 
reionization, it is possible that a large fraction of pre-reionization 'first-light' 
sources are insufficiently dusty and metal rich to be detectable as continuum 
sources. Low-metallicity galaxies should be detectable by fine-structure C and 
O far-IR line emission, however. It would be tremendously exciting to see the 
birth of the first metal-enriched dusty systems with ALMA, and so perhaps 
to determine directly the redshift limit for submm surveys. Of course, even if 
this were possible, ALMA would still have a long and fruitful career studying 
the detailed astrophysics of galaxies out to and beyond redshift 5, while the 
search for the most primitive galaxies in the second and third decades of the 
century is taken up by space-based mid-IR interferometers and the SKA radio 
telescope (see Fig. 8). 

Submm observations of the distant Universe are a new tool for probing the 
earliest and most dramatic stages of the evolution of galaxies. Over the years 
to come, the capabilities of submm-wave observatories, and our understanding 
of the Universe in this new window, should continue to advance dramatically. 
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